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ABSTRACT 

The Oort constants describe the local spatial variations of the stellar streaming field. The classic 
way for their determination employs their effect on stellar proper motions. We discuss various problems 
arising in this procedure. A large, hitherto apparently overlooked, source of potential systematic error 
arises from longitudinal variations of the mean stellar parallax, caused by intrinsic density inhomogeneities 
and inter-stellar extinction. Together with the reflex of the solar motion these variations by mode mixing 
create contributions to the longitudinal proper motions Jl^ (£) that are indistinguishable from the Oort 
constants at < 20% of their amplitude. Fortunately, we can correct for this mode mixing using the 
latitudinal proper motions IJ>i(£). 

We use about 10 6 stars from the ACT/Tycho-2 catalogs brighter than V ~ 11 with median proper 
motion error of w 3masyr . Taking every precaution to avoid or correct for the various sources of 
systematic error, significant deviations from expectations based on a smooth axisymmetric equilibrium, 
in particular non-zero C for old red giant stars. We also find variations of the Oort constants with the 
mean color, which correlate nicely with the asymmetric drift of the sub-sample considered. Also these 
correlations are different in nature than those expected for an axisymmetric Galaxy. 

The most reliable tracers for the "true" Oort constants are red giants, which are old enough to be 
in equilibrium and distant enough to be unaffected by possible local anomalies. For these stars we find, 
after correction for mode-mixing and the asymmetric drift effects, A ss 16, B « —17, A — B ss 33, and 
C w — lOkms -1 kpc _1 with internal errors of about 1-2 and external error of perhaps the same order. 
These values are consistent with our knowledge of the Milky Way (flat rotation curve and £1 = A — B as 
28 ± 2). Based on observations made with the ESA Hipparcos astromctry satellite. 

Subject headings: Stars: kinematics - Galaxy: fundamental parameters - Galaxy: kinematics and dynamics - 
Galaxy: solar neighborhood - Galaxy: structure 



1. Introduction 

For over eight decades, the distribution and kinemat- 
ics of stars in the solar neighborhood has been studied 
in order to infer the structure of the Milky Way galaxy. 
Kapteyn & van Rhijn (1920) used star counts to de- 
termine the size and thickness of the Milky Way. Fur- 
thermore, by assuming hydrostatic equilibrium perpen- 
dicular to the Galactic plane, the radial velocities and 
proper motions of nearby stars allowed Kapteyn (1922) 
to make the first reasonable estimate of the mass of the 
Milky Way. However, Oort (1927a) pointed out that the 
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mass of Kapteyn's Galaxy is not large enough to keep 
the globular clusters and RR Lyrae stars bound to the 
Galaxy. Lindblad (1927) proposed that the sub-systems 
of high-velocity stars and globular clusters as well as 
that of nearby low-velocity stars have the same axis of 
symmetry and a common center - that of the globular 
clusters (Shapley 1918). Lindblad further hypothesized 
that each of these sub-systems is in dynamical equilib- 
rium, and that the sub-systems with the largest amount 
of rotation will have the flattest space distribution and 
the smallest peculiar velocities (Jeans 1922), in agree- 
ment with the then available data. 

The motions of the stars in the solar neighborhood 
can be interpreted as a streaming (average) velocity plus 
random motions. In disk galaxies, the first dominates 



over the second: such stellar systems are said to be dy- 
namically cold. In the solar neighborhood, for instance, 
the velocity dispersion in the plane, i.e. the amplitude 
of random motions, is 45kms~ 1 for the old stellar disk 
and 18kms _1 for early-type stars, while the streaming 
motion is of the order of 200 km s -1 . In the cold limit 
of vanishing random motions, the streaming is along the 
closed orbits supported by the gravitational potential. 
Thus, knowledge about the potential of the Milky Way, 
and hence its mass distribution, can be gained from stud- 
ies of the stellar streaming velocities. 

Oort laid the theoretical basis of this method with his 
pioneering paper (1927a). To begin with, let us follow 
Oort and consider the cold limit in which all stars move 
on closed orbits. Oort himself actually considered the 
Milky way to be axisymmetric, but his analysis is eas- 
ily generalized to non-axisymmetry (Ogorodnikov 1932; 
Milne 1935; Chandrasekhar 1942). At each position x in 
the Galaxy, there is a unique streaming velocity v (here, 
we ignore the possibility of orbit crossing). The differ- 
ence Sv between the velocity at some point in the Galaxy 
and that at the Sun may be expanded in a Taylor series 
(with local Cartesian coordinates: e x and e y pointing in 
directions £ = 0° and £ = 90°) 



Sv = H • x + 0(x 2 ) 



(1) 



with 
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The parameters A, B, C, and K are the Oort constants, 
they measure the local divergence (K), vorticity (£?), 
azimuthal (A) and radial (C) shear of the velocity field 
generated by closed orbits. A star at Galactic longitude 
£ and distance d from an observer and moving with ve- 
locity v relative to the latter, is observed to have radial 
and tangential velocity 



Vd = v x cos £ + v y sin£ 
vg = —v x sin £ + v y cos £. 



(3a) 
(3b) 



Inserting equation (1) with x = d(cos£,sin£), whereby 
assuming the observer participates in the streaming field, 
one finds 

Vdd- 1 = K + Asin2£ + C cos2£ + 0([d/ R ] 2 ) (4a) 
Vid- 1 = B -Csin2£ + Acos2£ + O([d/R ] 2 ).('ib) 

It is worth emphasizing that terms quadratic in d/Ro 
contribute to the m — ±1,3 harmonics e lml in equa- 
tion (4) and thus do not interfere with the Oort constants 
- however, terms of higher odd orders do. Similarly, the 
deviations of the solar velocity from the local streaming 
velocity leads tom = ±l harmonics, see below. 

The Oort constants may also be expressed in terms 
of cylindrical coordinates (R, tp) with the Sun at (Rq, 0) 



(cf. Chandrasekhar 1942) 
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(5a) 
(5b) 
(5c) 
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evaluated at the solar position (we use the convention 
that tp increases clockwise, i.e. in the direction of Galac- 
tic rotation). In the axisymmetric limit, C = K = 0, 
and 
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equivalent to the equations actually given by Oort. 
In an axisymmetric Galaxy, the circular (closed) orbits 
have velocity v 2 = R d$/dR, and measurements of the 
Oort constants provide a direct constraint on the Galac- 
tic potential $. For instance, a harmonic potential re- 
sults in solid-body rotation, A = 0, and B equal to the 
(constant) rotation frequency, a flat rotation curve gives 
A = —B, and the case of all the mass concentrated in the 
inner Galaxy yields A = —3B. From the then available 
radial velocities and proper motions, Oort (1927b) found 
A w 19kms _1 kpc _1 and B w -24 kms -1 kpc -1 (with 
large uncertainties, though). This was clear evidence 
for a rotation of the Milky Way (not well established at 
that time), and ruled out Lindblad's suggestion that the 
Milky Way rotates like a solid body. Note, however, that 
zero C and K do not necessarily imply axisymmetry of 
the Galactic potential; alternatively, the Sun might be 
located near a principal axis of an elliptic potential (Kui- 
jken & Trcmainc 1994). 

Given their pre-eminent importance, it should not 
come as a surprise that the observational determination 
of the Oort constants has been high on the astronomers' 
priority list. However, measuring these streaming veloci- 
ties directly is no simple task, mainly because of the lack 
of an appropriate reference system which is unaffected 
by systematic motions and not involved in the Galac- 
tic streaming. For example, stellar positions and proper 
motions in the fundamental stellar catalogs (e.g. AGK3, 
FK4, FK5) are based on transit time measurements. The 
so-determined proper motions are absolute, but with re- 
spect to a reference frame tied to the Earth and the 
orbits of Solar system objects, and are thus useful to 
determine Earth's precession rate and the motion of 
the equinox (e.g., Fricke 1971; Clube 1972; Miyamoto & 
Soma 1993; Lindgren & Kovalevsky 1995). In order de- 
termine the Galactic rotation rate, it is more useful to de- 



2 



termine proper motions relative to an extra-galactic iner- 
tial reference frame. In more modern proper motion pro- 
grams, distant galaxies have been used for this purpose 
(e.g., NPM, SPM, KSZ, Bonn, Potsdam, see Kovalevsky 
et al. 1997, for a summary). However, until recently their 
usefulness has been limited by the so-called "magnitude 
equation" problem 3 . Comparison with the Hipparcos 
data showed that several of these extra-galactic reference 
frames were not exactly inertial. Some significant resid- 
ual spins (0.25 - 1.25 mas yr _1 ) of the pre-Hipparcos co- 
ordinate systems were found (see Kovalevsky et al. 1997, 
for a review) , which can lead to systematic errors in the 
Oort constant B of 1.2-6kms _1 kpc -1 . 

So far, we have only considered the cold limit and ig- 
nored the effect of random motions. However, we will 
see that this is not suitable for most stars in the so- 
lar neighborhood. That is, the mean velocity field of a 
group of stars differs systematically and significantly in 
its divergence, vorticity and shear from the (hypothet- 
ical) velocity field induced by closed orbits. Moreover, 
higher-order terms in the Taylor expansion and other ef- 
fects imply that the Fourier coefficients in equation (4) 
are not identical to the divergence, vorticity, and shear 
of the mean velocity field, let alone the closed-orbit ve- 
locity field. These problems (discussed in some detail 
in §2 with particular emphasize on proper-motion data) 
introduce significant systematic errors in the values de- 
rived for the Oort constants. In most previous studies, 
many of these problems have been ignored, a fact that 
may well explain the diversity among the values derived 
for the Oort constants from different data and by dif- 
ferent methods (see Kerr & Lynden-Bell 1986; Oiling & 
Mcrrifield 1998, for reviews). The main objective of the 
present paper is, after having recognized these system- 
atic errors, to avoid them as far as possible both by a 
careful analysis of the data and a careful interpretation 
of the results. §3 details our analysis technique. In §4, 
we motivate our selection of the ACT catalog and ana- 
lyze its proper motion data. The results are discussed 
and interpreted in §5. Finally, §6 sums up and concludes. 

For the convenience of authors and readers, the units 
used throughout this paper will be consistent and will 
not always be explicitly given. Distances are measured in 
kpc, velocities in kms -1 , and frequencies, such as proper 
motions and the Oort constants, in km s -1 kpc" 1 . Stel- 
lar parallaxes, denoted by the symbol ir, are considered 
to be inverse distances and, consequently, have dimen- 
sion kpc -1 , which corresponds to measuring the actual 
parallax angle in milli-arcseconds. Luminosities and col- 
ors are given in magnitudes as usual. 



3 Essentially, accurate position determination is hampered by the 
fact that stellar images are non-spherical as a result of telescope 
tracking errors and the non-linearity of the photographic plates 
(van Altena, private communications) 



2. The Oort Constants in Practice 
2.1. Deviations from the Cold Limit 

As already pointed out by Oort (1928) and later dis- 
cussed by Kuijken & Tremaine (1994), only in the cold 
limit of vanishing random motions can we interpret the 
mean streaming velocity v as the velocity v of closed 
orbits supported by the Galactic potential. In general, 
there is a systematic difference, and we may write 



v - v„ 



(7) 



with the asymmetric drift velocity v a . The asymmetric 
drift expresses the lag of the mean velocity with respect 
to the local closed orbit. To explicitly make this distinc- 
tion, let us write in analogy to (7) 



A — A A a , 
B = B — B a , 
C = C-C a , 
K = K-K n . 



(8a) 
(8b) 
(8c) 
(8d) 



Here, A, B, C, and K represent the mean velocity field 
(its divergence, vorticity and shear) of a group of stars 
and might be evaluated from equation (2) or (5) by re- 
placing v with v, while A a , B a etc. follow upon replacing 
v with v a . 

For an axisymmetric Galaxy, there is only an az- 
imuthal component v av , which, for random motions 
much smaller than the rotation velocity, can be approx- 
imated by Stromberg's relation (see Binney & Tremaine 
1987, eq. 4.34, for a derivation from the Jeans equations) 



2v v v av ~ a R 



_ 1 _ 



dlnR 



■ (9) 



Here, <r 2 is the (square of the) velocity dispersion tensor 
and p the stellar density. Thus, the asymmetric drift 
is a function of the radial velocity dispersion but also 
depends on the circular velocity, the axis-ratio of the ve- 
locity dispersion ellipsoid, as well as on gradients in the 
dispersion and stellar density. We might use this rela- 
tion to estimate the expected effect on Oort's A and B. 
We refer the reader to Lewis & Freeman (1989) for a 
determination of the radial dependence of the radial and 
tangential velocity dispersions. We proceed by neglect- 
ing the radial variation of the last term in (9) and by 
assuming that p and er 2 vary exponentially with scale 
lengths Rd and R a 2 . We then find: 



d\nv aip ^ 1 k ( 1 1 
dR ~ ~R^ + 2^ R^ 



d In? 



dR 



(10) 



where we have used k v aip ~ a R with k — (80 ± 5) kms 1 
(Dehnen & Binney 1998, hereafter DB98). 

For R a 2 — 0.45i?o (Lewis & Freeman 1989), one finds 
A a rts 0.14v aip and B a Q.Qlv aipi almost independent 
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of Rd, Ro, and v v . That is, B is hardly affected, while 
A a might be as large as 3 for red stars (using v atp of 20 
according to DB98's findings). 

This derivation of A a and B a is based on the assump- 
tion of dynamical equilibrium, and hence not appropri- 
ate for young stellar populations, whose lumpy phase- 
space structure (moving groups and inhomogeneous spa- 
tial distribution) indicates non-equilibrium. Further, 
note that the numerical values derived above from equa- 
tion (10) are only valid in the Solar neighborhood as k 
and v av were determined locally. 

It is important to notice that A, B, C, and K are func- 
tions of Galactic height z (already because the asymmet- 
ric drift depends on z). One might think that one could 
just extend the Taylor expansion (1) into the vertical 
direction. However, this does not work, mainly because 
the scale height of the stellar disk, and hence any pos- 
sible kinematic gradients, is much smaller than reason- 
able sample depths, i.e. high-order terms are significant. 
(Also, for symmetry reasons, the first non-trivial term is 
quadratic and its distance dependence would not nicely 
vanish in the resulting proper motions.) 

2.2. Deviations from Axisymmetry 

For an axisymmetric Galaxy, as originally considered 
by Oort, the closed orbits are circular and at each ra- 
dius there exists exactly one such orbit. However, the 
Galaxy is not axisymmetric, it has a central bar and its 
disk appears to have 4, or a 4+2, armed spiral pattern 
(Vallee 1995; Amaral & Lepine 1997), but the infrared 
light, originating mainly from old stars, is dominated by 
a 2-arm mode (Drimmel 2000; Drimmel & Spergel 2001). 
As a consequence, the closed orbits are no longer circu- 
lar but elliptical. The orientation of these elliptic orbits 
changes when one crosses the co-rotation or inner and 
outer Lindblad resonances (CR, ILR, & OLR) of the bar 
or spiral pattern. Certain places near these resonances 
are visited by two or more differently oriented closed or- 
bits. It has even been proposed that the Sun is precisely 
at such a position on the OLR (Kalnajs 1991). How- 
ever, it seems now more likely from detailed models of 
the gaseous and stellar motions in the inner parts of the 
Milky Way that we are about 1 kpc or even less outside 
the OLR (Englmaicr & Gerhard 1999; Fux 1999; Dehnen 
1999, 2000). The spiral pattern of the Milky Way also 
imposes non-axisymmetric perturbations on the veloc- 
ity field (Lin, Yuan & Shu 1969; Mishurov Pavlovskaia 
& Suchkov 1979). Analyses of radial velocities and now 
also of proper motions of several young tracer popula- 
tions indicate that the Sun is located close to the radius 
of co-rotations of the spiral density wave (Ai? = ±1 
kpc, AQ — ±2kms~ 1 kpc -1 ; Creze & Mennessier 1973; 
Amaral & Lepine 1997; Mishurov & Zenina 1999). 

Clearly, the Oort constants, defined in terms of the 
closed-orbit velocity field, are ill-defined at these posi- 
tions, and will behave discontinuously when crossing one 



of the three major resonances. Of course, in reality one 
measures the Oort constants from stars that are not on 
closed orbits. In this case, the discontinuities of the Oort 
constants will be replaced by a more gradual transition. 

2.3. Measuring the Oort Constants 

The Oort constants are measured by determining the 
Fourier coefficients in equation (4). However, in real- 
ity the Galaxy is not flat but three-dimensional, and we 
cannot use this equation directly. To generalize for stars 
out of the plane, we use 

n = d~ 1 cos b 
v r = Vd cos b 
~p t * =irve 
fi b = -7T v d sin b, 

where n denotes the star's parallax, d the distance pro- 
jected onto the Galactic plane and 7% = Hi cos& (the 
quantity actually measured in astromctry). Moreover, 
the Sun is not moving with the local streaming velocity 
but with some velocity Vq = (U, V, W) relative to it. 
We thus finally have for the observable proper motions, 
neglecting higher-order terms in equation (1), 

7% = ~pu sin £ — ~p v cos £ 

+ cos b (B — C sin2£ + A cos 2£) (11a) 
~p b = sin b (jLy cos £ + ~p v sin £) — ~p w cos b 

- sin b cos b (K + A sin 2£ + C cos 2£) , (lib) 

with the reflex of the solar motion 

7*© = (~Pu> ~Pv >7V) = kvq. (12) 

Note that we have written A instead of A etc.; the sym- 
bols with a tilde are defined by equation (11), i.e. as 
Fourier coefficients of the mean proper motion for stars 
at the same distance, while those with a bar are the di- 
vergence, vorticity, and shear of the mean velocity field 
for a group of stars. Strictly speaking, these are two dif- 
ferent sets of quantities, and we can only hope that they 
are not too different. Of course, all we can hope to mea- 
sure are A, B, C, and K, while we want A, B, C, and 
K (actually, A, B, C, and K). We now discuss potential 
sources of discrepancies between these two sets of quan- 
tities, and further problems in measuring the Fourier 
coefficients. 

2.4. Systematics with Sample Depth 

Equations (5) are based on the first-order Taylor series 
(1) and neglect higher-order terms, which might become 
important for deep samples. Including the next order 
adds six new unknowns, while only the m —1, 3 Fourier 
terms in 71^* are affected, i.e. only two more constraints 
are available (the cos 3^? and sm3£ terms). Thus, the 
next higher-order expansion is not generally soluble from 
proper motion data alone. 
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For illustration, we now consider the effect of a purely 
axisymmetric model, which adds only one additional un- 
known per Taylor order (Pont, Mayor & Burki 1994; 
Feast & Whitelock 1997). To next order in d/Ro, we 
then get at 6 = 0: 



(M**)o 

d 
8Rq 



d 
8i?o 

d 2 v, 

-»0 



6A 



OR 2 



OR 2 

6 A ) cos3£, 



where (jx^*)o refers to the expression (11a). From anal- 
yses of Cepheid kinematics, the value of d 2 v ip /dR 2 \R is 
estimated to be smaller than 3kms _1 kpc -2 in modu- 
lus. With this estimate and A w 15, one finds that the 
term 0(x) contributes up to —2.5 and 1.8 to the cosl 
and cos3£ Fourier terms, respectively, for d/Ro = 1/8. 
Thus, even for a rather smooth streaming field higher- 
order terms may become important already at modest 
sample depth. 

2.5. Effects of a Non-Smooth Streaming Field 

Because of local anomalies in the Galactic force field, 
e.g. caused by spiral arms and other distortions, the 
streaming field inevitably exhibits small-scale oscillations 
on top of an underlying smooth field. These oscilla- 
tions give rise to significant higher-order terms in the 
Taylor expansion (1). Since in practice the Oort con- 
stants are measured from the kinematics of large stellar 
samples with a finite depth, these high-order terms be- 
come important. In the present context, however, we are 
predominantly interested not in the small-scale but the 
large-scale behavior of the streaming field. When using 
a deep sample, i.e. a big volume, one might hope that 
the small-scale wiggles average out and one is left with 
the Oort constants due to the smooth equilibrium field. 
We have simulated this effect in Figure 1, where an ax- 
isymmetric streaming velocity (black line in upper panel) 
with wiggles of wave-length A w 2 kpc and amplitude of 
only two percent of an otherwise smooth (power-law) ro- 
tation velocity is assumed. From this model, and the 
exact proper motion equation 



fii(£,d) 



v tp (R) [Rq cos I — d]/R — Vip(Ro) cos I 



(13) 



we measured the Fourier coefficients A and B as function 
of distance d (Fig. lb). Assuming that the A and B 
obtained in this way measure the local streaming field, 
one may approximate that locally 

v v ~ 2R A- R(A + B). 

The straight lines in Fig. la, reaching from Rq — d to 
Rq + d, represent these approximations. One finds that 
for nearby stars (d <C A) the measured Oort constants A 
and B accurately represent the (wiggly) local streaming 
field. For distances d similar to or larger than A, A and B 
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Fig. 1. — Effect of a wiggly streaming velocity on the measured 
Oort's constants. The curved (black) line in (a) represents an ax- 
isymmetric model for the streaming (rotation) velocity. From the 
exact d) predicted by this model, we have measured, as func- 
tion of distance, A and B as Fourier coefficients (b). Recovering 
the rotation velocity from A and B obtained in this way yields the 
straight lines in (a). The radial extent of these lines corresponds 
to twice the distance. 



represent the underlying smooth part. For intermediate 
distances, there is a transition with a range in Rq(A — 
B) [the estimate for v v (Ro)] comparable to twice the 
amplitude of the wiggles. The estimated gradient of v v 
(the slope of the straight [red] lines) is a strong function 
of d, giving sensible results only for d > A. 

The importance of these effects is likely to be 
strongest for young stars and weakest for old stars, be- 
cause their larger velocity dispersion makes old stellar 
populations less susceptible to small-scale features in 
the force field. From the above it is clear that even the 
smallest degree of non-smoothness (2%) in the Galactic 
force field can have a significant effect on the appar- 
ent values of the Oort constants (30%). In fact, one 
can show that the degree to which a stellar popula- 
tion with radial dispersion ctr responds to small-scale 
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perturbations is inversely proportional to a\ (Mayor 
(1974)). Thus, while young stars are ideally suited to 
trace Galactic fine-structure, old stars with their larger 
dispersion are much less influenced by any wiggles in 
the force field and are thus more suitable to study the 
large-scale potential of the Milky Way. A quantitative 
estimate for the critical wave-length to which a group 
of stars with dispersion or is just sensitive, may be 
given by the average epicycle diameter \/&<7r/k. With 
a epicycle frequency of k rs 36kms~ 1 kpc~ 1 and o~r of 
w lSkms -1 and rts 38kms~ 1 for blue and red stars 
(DB98) , we find a limiting wave-length of ~ 1 kpc and 
w 3 kpc, respectively. From the structure in the inter- 
stellar medium (ISM), we expect wiggles on the scale of 
about 2 kpc (Oiling & Merrifield 1998), which can result 
in 30% "errors" in the Oort constants. 

Groups of stars with a radial velocity dispersion larger 
than about 25kms _1 have an epicyclic diameter larger 
than the wavelength of the ISM-induced wiggle of the 
rotation curve. Main-sequence stars bluer than B — V = 
0.42 have dispersions smaller than 25kms _1 (DB98) and 
main-sequence lifetimes less than about 2 Gyr. Thus, 
the kinematics of stars younger than 2 Gyr may be sig- 
nificantly influenced by low-amplitude, small-scale per- 
turbations in the Galactic potential. Non-axisymmctric 
perturbations -which were not considered by OM98- 
could lead to even larger differences between the appar- 
ent Oort constants and the "true" Oort constants. 

In fact, the situation depicted in Figure 1 is likely to 
be a simplification. In the calculations leading to Fig. 1 
we assumed that the Sun is located half-way between the 
extrema of the rotation curve wiggle. That is to say, we 
assumed that the Sun is located in a special place with 
respect to the inevitable small-scale oscillations of the 
Galactic force field. Oiling & Merrificld's (1998) work 
suggests that, indeed, the Sun is not located in a "sweet 
spot" of the wiggly v v (R) curve (see also Amaral et al. 
1996). In that case, the widely employed technique of 
expanding the rotation curve and the velocity field to 
low-order (1) is not adequate. 

And finally, since there is evidence that the Milky 
Way exhibits significant azimuthal asymmetry (§2.2), 
the arguments presented above are likely to be over- 
simplifications. Nevertheless, the reasoning above gives 
us an indication as to the level of possible systematic 
errors in our analyses (see also §2.8.2). 

2.6. Correlations of Parallax with Kinematics 

In reality, of course, the stars are not just at one dis- 
tance, but distributed over the sampling volume. That 
is, for each £, one averages along the line of sight. In 
order to arrive at an equivalent to equation (11) for a 
sample of stars are various ir, one must assume that 
ttU = WU and analogously for ttV and nW, i.e. that 
parallaxes and velocities are uncorrelated. However, for 
stars associated with spiral arms one expects such cor- 



relations. A systematic error is introduced when using 
this assumption for stellar samples for which it is not 
justified. 

2.7. Mode Mixing 

When measuring the Oort constants from proper mo- 
tions, one often uses a magnitude-limited sample with 
little information about the stellar parallaxes. One then 
has to replace tt in (12) by the line-of-sight averaged 
parallax W (or, more precisely, ttvq by nv Q ). Due to 
non-uniformity both of the stellar space-density and of 
the extinction, W will inevitably depend on longitude. We 
might expand it into a Fourier series in I 

W = 7f ^1 + 2 c mn cos ml + Snm sin m£j , (14) 

where the pre-factor of two in front of the Fourier sum 
has been introduced for later convenience. Inserting (14) 
and (12) into equation (11) gives 

M£* = ( Bcosb + Snlji m - C^ijiyo) 

+ cost ( - 7iy + s^V-uo - c^/Jyo) 
+ sin^ ( ~p m - CniTi m - Stt27Vo) 
+ cos2£ ( A cos b - Sn-JijjQ - Ck+Tlvq) 
+ sin 21 ( - (7cos6 + c^-TI^o - s^+Jty^) 

+ 

(15a) 

~p b = ( - cos b [JI W0 + sin bK\ 

+ sin 6 [c^iUm + s n ip vo ]) 
+ cos£ ( — 2c 7r i cos b~p wo 

+ sin b [JI U0 + c 7r2 7i; 70 + s^MyoD 
+ sin I? ( - 25^1 cos 671^0 

+ sin 6 [JI V0 + s^uq - c^Mvol) 
+ cos2£ ( - 2c n 2 cos b~p WQ 

+ sin b [- cos b C + c^+JIuq - s^-TVo] ) 
+ sin 2£ ( - 2s k2 cos b~p wo 

+ sin b [— cos b A + s^+Tiyo + c^^Ji V0 ] ) 

+ 

(15b) 

where c w± = c n i ± c n3 , s v ± = s w i ± s n3 , s^ 2 ±4 = 
s n 2 i Stt4 and 



(16) 



Thus, the m = and m = 2 Fourier coefficients of Jif* 
measurable for a stellar sample are not identical to A, 
B, and C, rather mode mixing leads to additional con- 
tributions from the solar reflex motion. Likewise, and 
contrary to the classical no-mode-mixing case, the So- 
lar reflex motion terms (m = 1) have contributions from 
7Iy andJI U0 . 

2.7.1. The Size of the Effect 

It is instructive to estimate the possible size of the 
effect under realistic conditions like those we will en- 
counter in our application to proper-motion catalogs be- 
low. Let us assume that we have a stellar sample at low 
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latitudes with mean parallax W of 2 mas, correspond- 
ing to a typical depth of about 500 pc. Then m = 1,3 
amplitudes of the mean parallax of only 10% causes a 
contribution of about 1-3 km s" 1 kpc^ 1 to the observable 
m = 0, 2 harmonic of ~p£* , i.e. the Oort constants (with 
U w lOkms -1 , V ps 20kms~ 1 ). This is larger than 
the uncertainty of the raw Fourier coefficients. Thus, 
mode mixing dominates the error of the Oort constants 
when determined from proper motions surveys, unless 
the m = 1,3 amplitudes of W are much smaller than 
10% 4 . Most amazingly, the corresponding literature is 
absolutely void of any remarks on this nasty effect. One 



Proper Motior Variation (ACT Catalog) 




Longitude [degrees] 



Fig. 2. — The longitudinal variation for a group of "red" and 
"blue" stars. We plot ~p£* and /Ltj, in the top and bottom panels, 
respectively. The jagged lines with and without error bars repre- 
sent the red and blue main-sequence stars, respectively. We also 
plot the Fourier fits as the thick smooth lines, and list the number 
of stars used, as well as the reduced x 2 values. The curves with 
the larger amplitudes correspond to the red stars. Only stars with 
sin(6)| < 0.1 are included. 

may show that the longitudinal dependence of the mean 
parallax induced just by a smooth exponential stellar 
disk is small 5 . However, as we will see below, the stars 
are very non-uniformly distributed in I. This is a clear 
indication of inhomogeneity in the spatial distribution 
of the sample stars. This inhomogeneity is caused both 
by an intrinsic dumpiness of the stellar density and by 
extinction blocking the view through the Galaxy in a 
highly inhomogeneous and anisotropic way. Thus, we 
expect ?f to be non-uniform as well, causing considerable 
mode mixing. 



4 Of course, if the sample is deeper, the effect is alleviated. However, 

then other unwanted effects appear, see §2.5. 
5 For a (volume-complete) sample with depth d ma x in an exponential 

disk with scale length R^, s^m = and, to first order in d ma x/fioi 



This leads to mode-mixing contributions of 0.25 to 0.6 for early 
to late stellar types, respectively, independent of Wo (at 6 = and 
with R d = 2.5 kpc, U ss 10, and V between 10 and 25, DB98). 



In Figure 2 we present graphically the longitudinal 
variation of ~p e * and fib for two groups of stars with dif- 
ferent colors. The "jagged" lines represents the data, 
the smooth lines the Fourier fits. We clearly see that 
the vertical proper motion exhibits azimuthal variation 
of about ± 35% (bottom panel), indicating a changing 
average parallax with longitude. The reduced x 2 values 
listed are computed with respect to a 5 th order Fourier-fit 
model and indicate that the fits are "reasonable." How- 
ever, because the distribution of the proper motions in 
each longitude bin of 3° is not quite normal, the xt val- 
ues are only indicative of the goodness of fit. Assuming 
a simpler model with constant Hh(l) lead to 60% larger 
xl values, proving that the observed variation is real. 
The mean levels of the vertical proper motion indicate 
that the red stars are about three times closer than the 
blue stars. 

2.7.2. A Cure from the Effect? 

In order to correct for this effect, one needs some unbi- 
ased estimate of the relative Fourier coefficients (c ff , s T ) m 
for the mean stellar parallax, defined in equation (14). 

Fortunately, for b — equation (15b) simplifies con- 
siderably to 

~p b | b=0 = ~~PW0 ( 1 + 2 X! Cjrm C ° S S rrm Sin mA . 

m 

(18) 

Thus, for low latitudes we can use the vertical stellar 
proper motions to measure, apart from ~Pwq, the relative 
Fourier coefficients (ev, s^m of w and use them to correct 
the Fourier coefficients obtained from ~p ti , (see §3.2 for 
details of our measurement technique). However, this 
works only, if the vertical solar reflex motion Wq, or 
equivalently, the mean vertical motion of the stars with 
respect to the Sun is constant with £. A stellar warp, for 
example, would invalidate this assumption. Fortunately, 
it seems from Dehnen's (1998) analysis of the local stellar 
velocity distribution that the stellar warp only weakly 
affects our analysis: only stars at R — R > 1 kpc or 
with high V-velocity are affected. 

If this assumption of constant Wo is doubtful, or if 
high-latitude stars are concerned, the only possible cure 
of the problem is an unbiased, but not necessarily very 
accurate, distance measure, which can be used to in- 
dependently estimate the Fourier coefficients of W. Us- 
ing the photometric parallaxes for this purpose does not 
work - we tried it - presumably since extinction af- 
fects the distributions of true and photometric parallax 
in fundamentally different ways, but perhaps also be- 
cause of systematic errors in the photometric parallax 
(cf. Malmquist bias). 

2.7.3. Mode Mixing and Distance Effects 

In §2.4 we have seen that for a conventional param- 
eterizations of the average Galactic streaming field, dis- 
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tance effects result in cos £ and cos 3£ amplitudes of 
order of several kms _1 kpc _1 . Such contributions are 
comparable to those expected from mode-mixing effects 
(§2.7.1). Thus, for samples that are sufficiently deep for 
these distance effects to be important one cannot dis- 
entangle the two contributions, and a correction of the 
mode mixing from proper motion data alone becomes 
impossible. We will return to this point in the discus- 
sion of the results obtained from proper-motion data. 

2.8. Measurement Accuracy 

2.8.1. Internal Uncertainties 

It is instructive to estimate the expected uncertainty 
for A, B, and C determined from the proper motions of 
a large stellar sample (K can only be sensibly measured 
from radial velocities). The observed proper motion for 
an individual star can be described by the mean plus 
a random component and an error term. The mean as 
function of (l,b) is given by equation (11) (neglecting 
mode mixing). The random component has two sources: 
random stellar motions leading to a random component 
proportional to parallax, and scatter in stellar parallaxes 
resulting in scatter in the terms arising from the solar 
reflex motion. 

For a stellar sample uniformly distributed inside 
a radius c? max , with mean motion with respect to 
the Sun of 20kms~ 1 , and a one-dimensional veloc- 
ity dispersion in the plane of 30kms _1 , the disper- 
sion of the random proper motion amounts to roughly 
50/d max kms _1 kpc -1 . Thus, in order to measure the 
Oort constants to an accuracy of e, we need of the order 
of N = 2500 x d m l x x (e/kms^kpc- 1 )- 2 stars. If one 
increases the accuracy by using ever deeper samples, one 
must be aware that, as discussed in §2.5 above, the mea- 
surable Oort constants are functions of sample depth. 
From this estimate, it is clear that the measurement 
uncertainties in the observed proper motions are almost 
always negligible in the error budget. 

2.8.2. Systematic Errors 

There is, however, another well known source of un- 
certainty: a residual systematic rotation of the astro- 
metric reference system used to measure the individual 
proper motions. The rotation frequency simulates a vor- 
ticity of the velocity field, and its z component is indis- 
tinguishable from B. For ground-based proper- motion 
surveys, one can only with great difficulty establish an 
inertial reference frame, and in many cases some resid- 
ual rotation remains (see Kovalevsky et al. 1997, for a 
recent review). For example, the rotation of the Hip- 
parcos reference frame is estimated to be no larger than 
0.25 masyr -1 = 1.18kms _1 kpc" 1 (Kovalevsky et al. 
1997). This implies that B cannot be accurately mea- 
sured, even when using the best astrometric reference 
system currently available! 



3. The Measurement Technique 

The classical way used to determine the Oort con- 
stants from proper motions fi = (jig*,^) is a least- 
squares fit of equation (11) to the observed proper 
motions. Hanson (1987), for instance, employed this 
method to determine just A and B from the NPM survey. 
This straightforward approach, however, cannot be rec- 
ommended. The main problem is that the actual func- 
tional form of the mean proper motion ~p(£, b) is possi- 
bly not well described by the fitting function (see §2). As 
with all parametric methods, such a mismatch inevitably 
leads to systematic errors in the values obtained for the 
parameters. 

As already discussed in §2.1 and §2.3 above, the mea- 
surable Oort constants and the relative mean proper mo- 
tions ~pu, ~pv, Hw are functions of stellar type and lat- 
itude. These variations introduce an a priori unknown 
dependence of p on latitude b. In his analysis Hanson 
accounted for a variation of p Q , but did not find a lat- 
itude dependence of A and B. Additional to the ex- 
pected variations with latitude, there will be variations 
of ~p with longitude as discussed in §2.4 and §2.5. The 
stars are always non-uniformly distributed in I such that 
the harmonics e lml are not orthogonal under averaging 
over £. This implies that, if there is (non-anticipated) 
power in higher-order harmonics, fitting only the low- 
order (m < 2) harmonics again introduces systematic 
errors in the parameter values obtained. 

In order to avoid these problems, we expand ~p in 
a Fourier series including high orders and restrict our- 
selves to a small range of latitudes. This is an essentially 
non-parametric approach in that it allows for a general 
functional form for ~fi(l) and avoids variations at high 
latitudes. To be more specific we determine the coeffi- 
cients ctm, S£m, Cb m , and s bm of 

m max 

~fit* = cm + ^ ct m cos m£ + st m sin ml , (19a) 

m>0 

Mb = c b0 + Cbm cos ml + Sbm sin ml (19b) 

m>0 

by a least-squares fit to the observed ~p e * and Hb (see 
also §3.3.1 below). 

3.1. Weighting and Exclusion of Outliers 

In order to minimize the errors, one might weight in- 
dividual stars with the inverse variance of their proper 
motion. This is essentially given by the stars parallax 
times the velocity dispersion, see §2.8 above. The latter 
is a function of color, and, since we will analyze the stars 
in color bins, the weights reduce to l0 _0 - 4m where m is 
the apparent magnitude in some passband. Such a tech- 
nique greatly increases the importance of the faintest and 
most distant stars in the sample, i.e. those objects that 
are most likely affected by extinction. Experimenting 
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with such a scheme, we found no significant reduction in 
the estimated uncertainties. 

Another method to reduce the errors is to exclude 
high-proper-motion stars. For example, K-cr-clipping ex- 
cludes stars whose jL t „ or /U& deviates from the mean at 
their longitude by more than k times the dispersion a. 
This approach is somewhat dangerous as it uses kine- 
matic information itself in a study aimed at kinematic 
quantities. 

The most conservative technique is to exclude just 
bright stars, since these are most likely very near and 
their proper motions are dominated by random veloci- 
ties. 

3.2. Correction for Mode Mixing 

Having obtained the raw Fourier coefficients (q, si, Cb, s&) 
we might want to correct them for the mode mixing de- 
scribed in §2.7, i.e. extract the Fourier coefficients that 
would have been measured without longitudinal varia- 
tions in the mean parallax. Equating (19) to (15) one 
gets a linear relation between the raw Fourier coefficients 
(c£, si, Cb, Sb) m and the unknowns 6 

(Mt/o> Mvo> Mwo> A B, C, ...), (20) 
where the dots stands for possible higher-order terms. 

3.2.1. At Low Latitudes 

As pointed out in §2.7, at low latitudes sin b ~ 0, i.e. 
the Galactic-rotation induced higher-order terms in ~p b 
vanish, and we can use ~p b to directly measure the mode- 
mixing coefficients a = (c,r,s T ) m , see Equations (14) 
and (18). In §4 below, we will determine the unknowns 
£ up to order m = 4, and need the mode-mixing coeffi- 
cients and hence the raw Fourier coefficients of Jt b up to 
order m = 5. Thus, we have 20 (non-linear) equations 
(9 from ~p e * and 11 from ~p b ) for 20 unknowns (10 from 
£ and 10 from a). 

3.2.2. At High Latitudes 

At high latitudes, we cannot estimate a. directly from 
the proper motions. However, if one is only interested in 
the low-order coefficients (up to m = 2 in equation (15), 
the situation is not too bad, as there are 10 equations 
for 6 unknowns in £ and another 6 in a. Thus, even 
though the equations have no unique solution, there is 
no complete freedom in the possible mode mixing. In 
order to get a unique answer, one may postulate that 
mode mixing is minimal, i.e. seek the smallest a. that 
solves the equations. 



'We will not try to solve for K separately for two reasons. First, 
the only distinction between jSvKO an d sin 6 A" in equation (15b) is 
due to the mode-mixing terms, which are usually small. Second, in 
almost all cases sin b K is much smaller than ~p-[yo an d neglecting 
it does not introduce a significant error. 



Note that in the analysis presented in §4 we will re- 
strict ourselves to low latitudes and use the procedure 
outlined in §3.2.1. 

3.3. Numerics and Error Analysis 

3.3.1. The Raw Fourier Coefficients 

Let us, for convenience, denote the vector of Fourier 
coefficients (c, s) m by c and that of the harmonic basis 
functions by y. Then, we determine c by minimizing the 
sample average ((/z — c • y) 2 ) giving 

c= (vtovr 1 ■ ivy)- (21) 

The variance V c is computed via error-propagation from 
the variance of (fiy) , which is given by the standard for- 
mula 

■wt (w ® tiy) - fay) ® (tiy) , , 
v <w> = a?31 ' V 2 > 

3.3.2. Correction for Mode Mixing 

The Cj, i — l,b, are linearly related to our desired 
quantity £, 

Ci = M 4 • i, (23) 

where the matrices Mj are given through equation (15) 
and depend linearly on the mode- mixing coefficients a. 
Let 

i=i,b 

be the error-weighted deviation from this relation. Solv- 
ing (23) is equivalent to minimizing (24) with respect to 
£ and a. The variance is then given by the £-part of 
twice the inverse of the Hessian matrix of x 2 . 

4. Analysis of the Proper Motion Catalogs 

The coordinate system defined by the Hipparcos Cat- 
alogue (ESA 1997) is currently the best realization of an 
inertial reference frame (Kovalcvsky et al. 1997). It is 
thus expedient to re-determine the Oort constants using 
proper motion data based on this reference frame. Un- 
fortunately, the Hipparcos Catalogue itself is not a kine- 
matically unbiased catalog, so that special care has to be 
taken when inferring kinematic properties from it (Bin- 
ney et al. 1997). A kinematically unbiased subsample 
extracted from the Hipparcos Catalogue by DB98 con- 
tains about 14000 stars with a maximum and average 
distance of about 100 and 80 parsec. Using our estimate 
from §2.8 for the accuracy of the Oort constants, we ex- 
pect for this sample an uncertainty of the Oort constants 
of about 4kms _1 kpc -1 , which is unacceptably large. 

In order to get more accurate estimates for the Oort 
constants, one needs deeper and/or larger samples. In 
the present paper we will employ the ACT Catalog (Ur- 
ban, Corbin & Wycoff 1998), which is essentially mag- 
nitude limited, and hence kinematically unbiased. The 
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ACT Catalog is a combination of the Astrographic Cat- 
alog (Eichhorn 1974; Corbin & Urban 1988; Urban et al. 
1998) with a median epoch near 1904, and the Tycho 
Catalogue (ESA 1997), which has a limiting magnitude 
of about V = 11.5. It contains, for nearly 10 6 stars, two 
color photometry 7 and accurate astrometric data (posi- 
tions and proper motions) defined in the Hipparcos ref- 
erence frame. 

We have also analyzed the Tycho-2 catalog (Hog et 
al. 2000a). This catalog entails a re-analysis of the data 
of Hipparcos' Tycho instrument, the astrographic cat- 
alog and many other astrometric datasets (Hog et al. 
2000b) . This catalog contains about twice as many stars 
as ACT catalog, where the majority of the "extra" stars 
are fainter than Vr ~ 10, while its astrometry is equiv- 
alent to that of the ACT catalog (Hog et al. 2000a). 
It might be tempting to use the larger Tycho-2 cata- 
log instead of the ACT. However, a twice larger sample 
means a roughly 40% larger limiting distance, so that the 
distance-dependent effects and higher-order terms may 
become important (see §2.4). Instead we decided to use 
the ACT catalog, and verify the results by a similar anal- 
ysis of the Tycho-2 data over the ACT magnitude range 
(8.5 < V T < 11.5). Since the ACT and Tycho-2 re- 
sults are entirely consistent with each other, we will use 
the term "ACT/Tycho-2" to describe our data in the 
remainder of this paper. 

Compared to the Hipparcos sub-sample, the sample- 
depth and the number of stars is increased by factors of 
about ten and 70, respectively. Estimating the expected 
accuracy for the Oort constants, one realizes that we 
may even divide this sample into subsets and still get 
reasonably accurate values for A, B, and C. 

4.1. Longitudinal Inhomogeneities 

As discussed in §2.3, spatial inhomogcneity and ex- 
tinction are expected to bias the Oort constants and to 
confuse them via mode mixing. We are able, though, to 
correct for the latter effect provided we restrict ourselves 
to low latitudes. For this reason, since the stellar density 
is highest there, and because we are mainly interested 
in the Oort constants in the Galactic plane, we restrict 
our analysis to stars with | sin 6| < 0.1 (|6| < 5.74°). 
After further excluding stars brighter than tub = 8.5 
(§3.1), we are left with 192 546 stars in the color range 
B-V e [-0.2,4]. 



7 We use the photometric system as defined in the ACT catalog, 
which in turn was taken from the Tycho Catalogue. That is to say, 
when we speak about apparent magnitudes we refer to "Tycho" 
magnitudes. These Tycho magnitudes have been used to calculate 
an approximate color in the Johnson system [(B — V)j]. Note that 
these color transformations depends somewhat on spectral type 
and luminosity class. Hence, by necessity, the ACT colors are only 
approximate, especially in the region where there is a transition 
between predominantly main-sequence stars to mainly giants (i.e., 
the sub-giant region). For details on the Tycho photometric sys- 
tem we refer to The Hipparcos and Tycho Catalogues (ESA 1997, 



Table 1 
The Color Bins 
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Note. — The boundaries of the (observed) B — V color; number 
N of stars fainter than Vt = 8.5; relative overabundance in n(i) 
(see §4.1); Average inverse parallax [in kpc, see equation (26b)]; 
our best estimate of the intrinsic color; the main-sequence lifetime 
tms ( m Gyr; bins 7a-10 comprise non-MS stars); and the percent- 
age fy of stars younger than 1.5 Gyr. For the main-sequence stars, 
we estimate fy from the MS lifetime and a constant SFR, for the 
post-MS star, we employ the Sn(£) method, see §5.1.1 for details. 

In order to get a feeling for the importance of the 
spatial inhomogeneities, we binned the stars into 60 bins 
in B — V and considered their longitudinal frequency, 
n(£). Figure 3 plots n(£) for the color bins 1 to 10 
(see Table 1) used in the analysis below. Most promi- 
nent in the bluest color bin (#1) are three peaks near 
I ~ 75°, 240°, and 285°, which presumably are caused 
by stars in the Sagittarius-Carina (£ = 75°, 285°) and 
the Orion-Cygnus (£ = 240°) spiral-arms. In order to 
quantify how important these peaks are, we measured 
the relative number Sn(£) of stars in three windows of 
30° width that were centered on the three peaks, and 
subtracted 1/4, the expectation value for a uniform dis- 
tribution. Figure 4 plots this measure versus the mean 
color for 60 narrow color bins (open), as well as the 19 
bins used in the analysis below (solid). The importance 
of the peaks is largest at the blue end, where the stars 
are very young and not mixed or settled into equilibrium. 
Moreover, these stars are very bright and can be seen out 
to a few kpc, such that patchy extinction can best con- 
tribute to the apparent non-uniformity. Along the main 
sequence (MS), where the stars get older, hence better 
mixed, and fainter, n{£) steadily becomes more uniform 
until B — V ~ 0.7. Then the peaks become important 
again near B — V w 1, where the red clump dominates 
the sample, while the giants at B — V > 1.2 are nearly as 
uniformly distributed as are the dwarfs near B — V ~ 0.7. 
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Fig. 3. — The distributions n(£) (in arbitrary units) of the stars in the color bins 1 to 10, see Table 1 for the limits in B — V. 

4.2. Data Analysis 

4-2.1. Binning and Analyzing the Data 

Before the analysis in terms of the Oort constants, 
we split the data it into ten color bins, labeled 1 to 10. 
These bins have been chosen to group together stars with 
similar n(£). Small bins had to be avoided, in particu- 
lar at red colors, to yield reasonable errors for the Oort 
constants - large errors would render differences between 
the results of adjacent bins insignificant. In addition to 
these ten distinct bins, we consider 9 bins of similar size 
(labeled la to 9a) whose stars are taken half from each 
of the nearest primary bins. The color limits of all the 
19 bins are given in Table 1, which also lists the numbers 
of stars fainter than Vr — 8.5 contained in each bin. 

To reduce the errors, we also analyzed the data after 
applying a modest K-cr-clipping with k = 5 and excluding 
not more than 5% of the stars in each of 20 longitudinal 
bins (the only usage ever of such bins). The raw Fourier 
coefficients for ~p f i , and ~p b and their errors were estimated 
up to TOmax = 32 as outlined in §3.3.1 and are tabulated 
up to m = 5 in Tables 2 and 3 for all 19 color bins. Note 
that if there were no mode mixing, A, B, and C could 
be read of Table 2 as C2, Co, and — S2- 

The mode mixing caused by longitudinal variations of 
7f has been corrected assuming the Fourier coefficients of 
JI b are due to this effect alone (cf. §3.2.1 and §3.3.2). The 
resulting mode-mixing corrected Fourier coefficients are 
given in Table 4 and are discussed in the next section. 

4-2.2. Secular Parallax and Asymmetric Drift 

One might estimate the mean parallax ttq and the 
azimuthal asymmetric drift v aip from the reflex of the 
solar motion and the assumption that there are no radial 
and vertical components to the asymmetric drift, i.e. 

v e = (U , V +v av , W ). (25) 



The red- most of the 60 narrow bins at B — V = 1.95 
is somewhat odd, as it has less stars in those peaks than 
even for uniform n(£). The likely reason is that these 
are stars subject to severe extinction, which made them 
appear much redder than they really are (there are no 
stars with intrinsic B — V > 2), and restricts them to re- 
gions where high extinction has diminished the numbers 
in the other color bins. 
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Fig. 4. — The relative over-abundance of stars, Sn(£), in three 
30° wide windows centered on the peaks that dominate n(£) for 
the bluest stars. The open symbols are for 60 color bins while the 
solid squares represent the 19 bins used in the analysis; see also 
Fig. 3 for n(£) of these ten bins. 
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Table 2 

Raw Fourier Coefficients for ~p, t * 



bin 


CO 


Cl 


si 


C2 


S2 


C3 


«3 


C4 


s 4 


C5 


S5 


1 


-11.8±1.1 


-13.9±1.7 


9.8±1.5 


12.4±1.5 


-0.2±1.7 


-0.1±1.4 


-5.4±1.8 


0.7±1.7 


0.0±1.6 


-0.1±1.7 


1.4±1.6 


la 


-12.0±0.6 


-14.6±1.0 


12.1±0.8 


11.8±0.8 


1.1±0.9 


0.4±0.7 


-4.8±1.0 


0.4±0.9 


1.4±0.9 


0.1±0.9 


1.3±0.9 


2 


-11.9±0.4 


-15.4±0.6 


13.3±0.5 


13.4±0.5 


1.7±0.6 


0.9±0.4 


-4.4±0.7 


0.1±0.5 


1.2±0.6 


-0.7±0.5 


0.8±0.6 


2a 


-11.9±0.3 


-17.2±0.5 


14.0±0.4 


14.2±0.4 


1.8±0.5 


0.7±0.4 


-3.4±0.5 


-0.1±0.5 


0.5±0.5 


0.0±0.5 


0.4±0.5 


3 


-11.9±0.4 


-17.9±0.6 


14.0±0.5 


14.7±0.5 


1.8±0.6 


-0.0±0.5 


-3.9±0.7 


0.2±0.5 


0.8±0.7 


0.6±0.6 


-0.0±0.6 


3a 


-11.8±0.4 


-19.3±0.6 


14.3±0.5 


15.6±0.5 


1.3±0.6 


-1.0±0.5 


-3.6±0.6 


0.8±0.6 


0.1±0.6 


1.0±0.6 


0.6±0.5 


4 


-12.0±0.4 


-21.9±0.6 


15.8±0.6 


14.8±0.6 


0.9±0.6 


-0.7±0.6 


-3.2±0.6 


0.6±0.6 


-0.5±0.6 


1.2±0.7 


0.8±0.6 


4a 


-12.1±0.5 


-25.6±0.8 


17.4±0.6 


14.7±0.7 


1.4±0.7 


-0.4±0.7 


-3.3±0.7 


0.9±0.6 


-0.0±0.8 


0.6±0.7 


-0.2±0.7 


5 


-12.4±0.6 


-28.8±0.9 


20.3±0.7 


15.7±0.8 


1.6±0.8 


-1.6±0.8 


-3.2±0.9 


1.5±0.8 


0.1±0.8 


0.4±0.9 


-0.8±0.8 


5a 


-11.7±0.6 


-36.2±0.9 


24.5±0.7 


18.2±0.9 


1.6±0.8 


-1.9±0.8 


-3.4±0.8 


-0.3±0.8 


-0.4±0.8 


-0.8±0.9 


-1.0±0.8 


6 


-10.7±0.7 


-43.9±1.0 


26.5±0.9 


18.8±1.0 


2.0±0.9 


-1.9±1.0 


-2.8±1.0 


-O.lil.O 


0.5±0.9 


-1.7±1.0 


-0.6±0.9 


6a 


-9.6±0.9 


-63.7±1.3 


30.9±1.1 


17.7±1.3 


3.0±1.2 


-3.7±1.2 


-1.8±1.3 


1.7±1.2 


0.5±1.3 


-0.4±1.3 


0.8±1.3 


7 


-9.6±1.0 


-75.7±1.5 


32.1±1.3 


16.4±1.5 


1.9±1.4 


-5.2±1.4 


-2.4±1.4 


1.3±1.4 


-1.1±1.5 


-0.5±1.4 


0.9±1.4 


7a 


-9.8±1.4 


-81.1±2.2 


28.9±1.7 


11.6±2.2 


-5.7±1.9 


-7.7±2.1 


0.8±2.0 


0.5±2.1 


-3.5±2.0 


-2.1±2.0 


-0.1±2.1 


8 


-9.8±1.4 


-60.2±2.2 


23.5±1.7 


13.7±1.9 


-9.4±2.1 


-2.6±1.8 


4.7±2.2 


-0.8±2.0 


-5.3±2.0 


1.8±2.0 


1.9±2.0 


8a 


-12.2±1.1 


-45.9±1.8 


20.1±1.4 


14.2±1.5 


-5.4±1.8 


0.2±1.4 


1.3±1.9 


0.1±1.6 


-3.8±1.7 


3.7±1.6 


2.1±1.8 


9 


-13.6±1.0 


-40.0±1.6 


15.3±1.4 


13.7±1.3 


-0.5±1.7 


-0.8±1.3 


-0.6±1.7 


1.6±1.6 


-0.8±1.4 


0.7±1.6 


O.lil.d 


9a 


-12.8±1.1 


-38.6±1.7 


13.9±1.5 


14.2±1.3 


-1.9±1.9 


0.0±1.5 


1.5±1.8 


1.0±1.8 


-2.0±1.4 


0.5±1.9 


1.1±1.4 


10 


-12.2±0.9 


-30.5±1.3 


13.5±1.2 


13.8±1.0 


-3. Oil. 5 


1.2±1.2 


1.7±1.4 


-0.9±1.4 


-1.4±1.3 


0.4±1.3 


1.7±1.3 



with K-tr-clipp 



1 


-13.3±0.6 


-15.3±0.9 


11.8±0.9 


13.1±0.7 


0.5±1.0 


la 


-12.5±0.4 


-15.2±0.6 


13.0±0.6 


13.2±0.5 


1.5±0.7 


2 


-11.9±0.3 


-15.8±0.4 


13.6±0.4 


13.5±0.4 


1.7±0.5 


2a 


-11.9±0.2 


-17.0±0.4 


13.8±0.4 


13.7±0.3 


1.7±0.4 


3 


-12.0±0.2 


-18.2±0.4 


13.8±0.3 


14.5±0.3 


1.5±0.4 


3a 


-11.9±0.2 


-19.5±0.3 


14.1±0.3 


14.8±0.4 


1.7±0.4 


4 


-11.8±0.3 


-20.9±0.4 


16.1±0.4 


14.6±0.4 


2.4±0.4 


4a 


-11.5±0.3 


-23.1±0.4 


18.1±0.4 


15.4±0.4 


2.5±0.4 


5 


-11.0±0.3 


-25.6±0.5 


20.3±0.4 


16.2±0.5 


2.7±0.5 


5a 


-10.9±0.4 


-31.9±0.5 


23.7±0.5 


18.0±0.6 


3.0±0.5 


6 


-10.7±0.5 


-38.7±0.7 


25.1±0.7 


18.6±0.7 


2.6±0.7 


6a 


-9.9±0.6 


-53.8±0.9 


28.3±0.9 


17.6±1.0 


3.2±0.9 


7 


-9.6±0.7 


-63.5±1.0 


30.4±1.0 


18.1±1.1 


1.5±1.0 


7a 


-lO.lil.O 


-66.2±1.5 


27.1±1.3 


13.2±1.5 


-7.4±1.5 


8 


-9.5±0.7 


-42.2±1.2 


19.7±0.9 


15.7±1.0 


-5.6±1.2 


8a 


-11.3±0.5 


-34.2±0.8 


16.3±0.7 


16.4±0.7 


-0.2±0.8 


9 


-12.7±0.5 


-33.0±0.7 


14.3±0.6 


15.4±0.7 


2.0±0.7 


9a 


-12.6±0.4 


-31.6±0.6 


13.9±0.6 


14.4±0.7 


1.0±0.7 


10 


-12.5±0.4 


-27.4±0.6 


13.0±0.6 


14.2±0.6 


-1.0±0.6 



0.3±0.8 -4.8±1.0 2.3±0.9 0.7±1.0 1.3±1.0 2.2±0.8 

0.4±0.5 -4.2±0.7 1.0±0.6 1.0±0.6 0.1±0.6 1.8±0.5 

1.2±0.4 -3.9±0.5 0.1±0.5 1.1±0.4 -0.1±0.5 0.9±0.4 

1.1±0.3 -3.4±0.4 0.2±0.4 0.9±0.4 0.3±0.4 0.6±0.4 

0.2±0.3 -3.4d=0.4 0.3±0.4 0.2±0.4 0.6±0.4 0.6±0.4 

-0.8±0.3 -3.1±0.4 0.2±0.4 -0.3±0.4 0.4±0.4 0.9±0.4 

-1.1±0.4 -2.5±0.4 0.1±0.4 -0.3±0.4 0.0±0.4 0.9±0.4 

-0.5±0.4 -2.7±0.4 1.0±0.4 0.0±0.4 -0.2±0.4 0.3±0.4 

-1.0±0.4 -2.7±0.5 0.8±0.5 -0.1±0.5 -0.1±0.5 0.1±0.5 

-1.7±0.5 -2.5±0.6 0.2±0.5 -0.1±0.6 -0.5±0.6 -0.1±0.5 

-1.9±0.7 -2.1±0.7 0.4±0.7 1.4±0.7 -0.6±0.7 -0.8±0.7 

-3.7±0.9 -2.2±0.9 1.3±0.9 0.6±0.9 -1.6±0.9 -0.0±0.9 

-3.5±1.0 -1.3±1.0 1.7±1.0 -0.8±1.1 -0.8±1.0 -0.5±1.0 

-3.9±1.4 1.7±1.5 0.5±1.4 -4.3±1.5 0.2±1.4 -2.3±1.5 

-0.3±1.0 3.8±1.2 -2.2±1.1 -3.8±1.1 2.6±1.1 l.lil.O 

0.2±0.7 1.2±0.9 -0.1±0.8 -1.1±0.8 1.3±0.8 0.8±0.8 

-0.5±0.7 0.6±0.7 0.9±0.7 0.4±0.7 -0.1±0.7 0.6±0.7 

-0.0±0.6 1.4±0.7 0.3±0.7 0.1±0.7 -0.9±0.7 1.2±0.7 

0.9±0.6 2.1±0.6 -1.2±0.6 0.3±0.6 -0.9±0.6 1.4±0.6 
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Table 3 

Raw Fourier Coefficients for ~p b 



bin 


CO 


Cl 


si 


C2 


S'2 


C3 


S3 


C4 


s 4 


C5 


S5 


1 


-11.9±1.7 


-1.3±2.7 


-1.4±2.0 


-0.0±1.5 


-2.9±3.0 


2.9±1.2 


-2.4±3.1 


1.0±2.4 


-0.7±2.4 


4.2±3.1 


-0.0±1.4 


la 


-11.3±1.0 


0.3±1.6 


-0.7±1.1 


-0.1±0.9 


-1.6±1.7 


2.9±0.7 


-2.7±1.8 


0.4±1.3 


0.6±1.4 


1.2±1.7 


-0.7±0.8 


2 


-10.0±0.3 


1.9±0.5 


0.2±0.4 


0.5±0.4 


-0.1±0.4 


1.9±0.4 


-0.8±0.5 


0.1±0.4 


1.2±0.5 


0.3±0.4 


0.3±0.4 


2a 


-10.3±0.2 


1.0±0.3 


-0.0±0.3 


0.4±0.3 


-0.7±0.4 


1.2±0.4 


-0.7±0.4 


0.2±0.4 


1.4±0.4 


0.5±0.3 


0.7±0.4 


3 


-11.0±0.4 


1.6±0.7 


-0.9±0.5 


-0.2±0.5 


-0.9±0.7 


1.1±0.4 


-1.4±0.7 


-0.0±0.4 


2.5±0.7 


0.5±0.6 


0.2±0.6 


3a 


-11.2±0.4 


2.3±0.6 


-0.8±0.5 


-1.0±0.5 


-0.8±0.6 


0.7±0.4 


-1.1±0.6 


-0.2±0.5 


2.5±0.6 


0.1±0.6 


0.6±0.5 


4 


-11.9±0.3 


1.9±0.5 


-0.0±0.4 


-0.4±0.5 


-1.7±0.5 


1.0±0.5 


-0.8±0.5 


-0.2±0.4 


1.4±0.5 


0.2±0.5 


0.5±0.5 


4a 


-13.4±0.3 


2.4±0.6 


0.2±0.4 


-0.0±0.6 


-1.7±0.5 


2.4±0.5 


-0.8±0.5 


0.0±0.5 


2.0±0.6 


0.6±0.5 


0.8±0.6 


5 


-15.2±0.4 


3.2±0.6 


-0.6±0.4 


-0.1±0.6 


-1.4±0.5 


2.7±0.6 


-0.6±0.5 


-0.0±0.5 


2.6±0.6 


0.0±0.5 


1.8±0.6 


5a 


-18.2±0.4 


3.7±0.6 


-0.4±0.4 


-0.8±0.6 


-1.5±0.5 


1.8±0.6 


-0.7±0.5 


-0.6±0.5 


1.7±0.5 


0.8±0.5 


1.4±0.5 


6 


-21.1±0.5 


3.0±0.8 


-0.5±0.5 


-2.2±0.8 


-1.9±0.6 


1.6±0.7 


-0.8±0.7 


-0.4±0.7 


1.0±0.7 


1.0±0.7 


1.0±0.7 


6a 


-25.9±0.6 


1.5±1.0 


-1.9±0.7 


-3.6±0.9 


-2.9±0.8 


0.6±0.8 


-2.2±0.9 


-0.1±0.8 


-0.4±0.9 


-0.2±0.8 


-0.7±0.9 


7 


-28.2±0.6 


1.2±1.1 


-2.3±0.8 


-4.7±1.0 


-3.7±0.9 


0.5±0.9 


-1.8±0.9 


-1.2±0.9 


-0.9±0.9 


-0.5±0.9 


0.2±0.9 


7a 


-28.4±1.0 


1.2±1.7 


-4.3±1.1 


-8.1±1.5 


-3.8±1.3 


0.0±1.5 


-1.6±1.3 


-1.7±1.4 


-1.1±1.3 


0.5±1.4 


2.8±1.4 


8 


-21.1±0.9 


2.3±1.6 


-5.2±1.0 


-5.6±1.3 


-1.3±1.3 


-0.8±1.3 


-3.9±1.4 


1.3±1.3 


-0.6±1.3 


1.9±1.3 


-0.2±1.3 


8a 


-15.4±0.7 


1.8±1.3 


-4.0±0.8 


-3.3±1.0 


-O.Oil.l 


-0.2±0.9 


-3.0±1.2 


-O.Oil.l 


l.lil.l 


1.5±1.1 


-0.7±1.0 


9 


-12.5±0.6 


0.5±0.9 


-3.9±0.7 


-2.9±0.7 


-1.1±0.9 


-0.3±0.8 


-2.9±0.9 


-1.0±0.8 


0.8±0.8 


2.5±0.8 


0.3±0.9 


9a 


-11.7±0.6 


0.9±1.0 


-4.1±0.8 


-2.7±0.7 


-l.Oil.l 


-0.4±0.8 


-3.0±1.0 


-0.3±1.0 


-0.1±0.9 


1.2±1.0 


-0.4±0.9 


10 


-9.5±0.6 


0.8±0.8 


-3.8±0.8 


-1.0±0.8 


0.6±0.8 


-0.8±0.7 


-3.5±0.9 


-0.4±0.9 


0.1±0.8 


-0.1±0.9 


-1.5±0.8 













with k 


-cr-clipping 












1 


-11.0±1.6 


-0.9±2.6 


-0.9±1.8 


-l.lil.l 


-3.2±2.8 


2.8±0.9 


-3.7±2.9 


1.3±2.2 


0.0±2.1 


2.2±2.8 


0.2±0.8 


la 


-9.8±0.3 


1.0±0.6 


0.4±0.4 


0.0±0.3 


-0.5±0.6 


1.9±0.3 


-1.0±0.6 


-0.1±0.5 


1.4±0.5 


0.0±0.6 


0.3±0.4 


2 


-9.9±0.2 


1.3±0.3 


0.1±0.3 


0.3±0.2 


-0.3±0.4 


1.6±0.2 


-0.4±0.4 


0.2±0.3 


1.1±0.3 


0.3±0.3 


0.4±0.3 


2a 


-10.2±0.1 


1.2±0.3 


-0.0±0.2 


0.4±0.2 


-0.3±0.3 


1.3±0.2 


-0.4±0.3 


0.3±0.2 


1.1±0.2 


0.5±0.3 


0.6±0.2 


3 


-10.5±0.1 


1.0±0.3 


-0.4±0.2 


0.4±0.2 


-0.8±0.2 


1.3±0.2 


-0.7±0.3 


-0.0±0.2 


1.3±0.2 


0.5±0.2 


0.5±0.2 


3a 


-ll.liO.l 


1.6±0.3 


-0.7±0.2 


-0.1±0.2 


-0.9±0.2 


1.3±0.2 


-1.0±0.2 


-0.2±0.2 


1.2±0.2 


0.2±0.2 


0.5±0.2 


4 


-12.3±0.1 


1.7±0.3 


-0.3±0.2 


-0.4±0.2 


-1.2±0.2 


1.5±0.2 


-1.2±0.3 


-0.1±0.2 


1.2±0.2 


-0.2±0.3 


0.5±0.2 


4a 


-13.3±0.1 


2.1±0.3 


-0.2±0.2 


-0.5±0.2 


-1.6±0.2 


2.1±0.2 


-1.2±0.3 


-0.1±0.2 


1.4±0.2 


0.1±0.2 


0.8±0.2 


5 


-14.7±0.2 


2.8±0.3 


-0.9±0.2 


-0.7±0.3 


-1.8±0.3 


2.0±0.3 


-1.2±0.3 


-0.1±0.3 


1.3±0.3 


0.6±0.3 


0.8±0.3 


5a 


-17.0±0.2 


3.7±0.3 


-0.4±0.3 


-0.6±0.3 


-2.0±0.3 


1.7±0.3 


-1.3±0.3 


-0.4±0.3 


0.8±0.3 


0.7±0.3 


0.6±0.3 


6 


-19.1±0.2 


3.7±0.5 


-0.2±0.3 


-1.2±0.4 


-2.6±0.4 


1.4±0.4 


-1.0±0.4 


-0.3±0.4 


0.6±0.4 


1.0±0.4 


0.6±0.4 


6a 


-23.8±0.3 


2.9±0.6 


-1.1±0.5 


-2.5±0.5 


-3.0±0.5 


0.5±0.5 


-2.1±0.5 


-0.2±0.5 


-0.4±0.5 


0.5±0.5 


0.4±0.5 


7 


-25.7±0.4 


3.8±0.7 


-1.7±0.5 


-3.5±0.6 


-3.5±0.6 


1.6±0.6 


-2.5±0.6 


-1.1±0.6 


-0.7±0.6 


1.2±0.6 


0.8±0.6 


7a 


-23.2±0.6 


1.8±1.0 


-4.7±0.7 


-5.2±0.9 


-2.9±0.8 


0.9±0.8 


-3.6±0.9 


-1.4±0.8 


-0.9±0.9 


2.5±0.8 


0.7±0.9 


8 


-15.1±0.4 


0.8±0.7 


-3.8±0.5 


-2.8±0.6 


-1.7±0.7 


-1.2±0.6 


-3.2±0.7 


0.2±0.6 


0.0±0.7 


1.2±0.7 


0.4±0.6 


8a 


-12.4±0.3 


1.7±0.5 


-2.1±0.4 


-2.0±0.4 


-1.2±0.5 


-0.1±0.4 


-2.2±0.5 


-0.6±0.5 


0.1±0.5 


0.6±0.5 


-0.3±0.4 


9 


-11.0±0.2 


1.7±0.4 


-2.4±0.3 


-1.6±0.4 


-1.3±0.4 


0.2±0.4 


-1.9±0.4 


-0.8±0.4 


0.3±0.4 


1.8±0.4 


0.1±0.4 


9a 


-9.6±0.2 


1.2±0.4 


-2.9±0.3 


-1.6±0.4 


-0.5±0.4 


-0.4±0.4 


-1.9±0.4 


-1.0±0.4 


0.1±0.4 


1.8±0.4 


-0.0±0.4 


10 


-8.4±0.2 


1.2±0.4 


-2.8±0.3 


-1.1±0.3 


0.8±0.4 


-0.5±0.3 


-3.0±0.4 


-0.7±0.4 


0.0±0.3 


0.8±0.4 


0.0±0.3 
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Table 4 

Results after Correcting for Mode Mixing 

























bin 




fJ'VO 




A 


B 


C 


C3 


S3 


c 4 


S4 



1 11.7±2.9 15.3±2.1 11.9±1.7 10.9±2.2 -11.6±1.5 -0.1±2.5 0.3±2.0 -3.6±2.7 -2.7±2.4 1.0±2.3 
la 13.4±1.7 15.5±1.3 11.3±1.0 8.4±1.4 -12.7±1.0 -2.0±1.5 1.6±1.2 -4.5±1.8 -1.1±1.3 4.8±1.6 

2 13.1±0.7 15.9±0.8 10.0±0.3 9.6±0.8 -13.3±0.6 -2.1±1.0 1.4±0.8 -5.1±1.0 -0.8±0.8 2.7±1.1 
2a 14.4±0.6 18.1±0.6 10.3±0.2 11.6±0.7 -12.9±0.5 -2.4±0.8 1.7±0.7 -3.9±0.8 -0.6±0.7 0.9±0.8 

3 15.0±0.7 18.3±1.0 11.0±0.4 12.0±1.0 -13.8±0.8 -4.0±1.4 2.6±0.7 -5.3±1.6 -O.Oil.O 2.2±1.5 
3a 15.7±0.8 19.0±0.9 11.2±0.4 12.7±1.0 -14.4±0.8 -4.1±1.2 2.4±0.9 -5.6±1.4 1.2±0.9 0.9±1.4 

4 17.8±0.9 22.7±0.8 11.9±0.3 11.3±1.1 -13.9±0.7 -2.4±1.0 2.2±1.0 -3.0±1.1 0.3±1.0 0.2±1.1 
4a 19.1±1.0 26.8±1.0 13.4±0.3 9.0±1.3 -14.4±0.8 -2.0±1.1 2.2±1.2 -3.7±1.2 -0.9±1.1 1.3±1.3 

5 21.8±1.0 29.7±1.1 15.2±0.4 9.8±1.4 -16.0±0.9 -3.2±1.1 1.5±1.3 -4.4±1.2 0.7±1.2 0.8±1.3 
5a 26.7±1.1 36.5±1.1 18.2±0.4 12.4±1.4 -15.9±0.9 -4.2±1.2 2. Oil. 3 -3.8±1.3 -1.3±1.3 -0.4±1.3 

6 30.2±1.3 43.1±1.3 21.1±0.5 13.8±1.8 -14.2±1.2 -4.4±1.5 3. Oil. 6 -3.2±1.6 -1.5±1.6 0.7±1.7 
6a 37.0±1.8 61.4±1.6 25.9±0.6 14.8±2.3 -12.8±1.5 -8.6±2.0 2.5±2.1 -0.2±2.1 2.3±2.0 4.7±2.3 

7 40.3±2.1 72.2±1.9 28.2±0.6 14.5±2.7 -12.9±1.8 -7.8±2.3 4.5±2.5 0.9±2.4 2.9±2.5 1.7±2.5 
7a 39.4±3.1 73.2±2.6 28.4±1.0 11.9±4.2 -14.4±2.7 -2.7±3.3 6.8±3.7 2.7±3.5 2.9±3.9 -5.4±3.7 

8 29.0±3.0 53.9±2.5 21.2±0.9 12.7±3.7 -16.4±2.5 -4.4±3.3 3.2±3.3 2.3±3.6 0.2±3.6 -2.1±3.6 
8a 22.8±2.5 41.4±2.1 15.4±0.7 12.7±3.1 -17.7±2.1 -5.5±3.0 5.7±2.6 -2.5±3.4 0.1±3.1 -0.0±3.1 

9 19.3±2.4 36.7±1.7 12.5±0.6 14.2±2.2 -17.4±1.7 -10.2±2.5 6.6±2.4 -1.6±2.6 0.9±2.5 0.6±2.4 
9a 17.6±2.8 35.3±1.9 11.7±0.6 14.2±2.3 -17.3±1.8 -9.8±3.0 5.4±3.0 1.6±3.1 1.8±3.1 2.0±2.7 

10 13.4±2.0 28.5±1.8 9.5±0.6 13.8±2.1 -16.2±1.5 -9.0±2.3 3.2±2.3 0.2±2.7 1.8±2.5 5.6±2.4 



with ft-cr-clipping 



1 15.1±3.0 16.7±1.7 11. Oil. 6 9.7±1.4 -13.3±1.2 -1.5±2.0 2.4±1.7 -4.1±2.0 1.1±1.6 3.8±1.7 
la 13.4±0.8 15.6±0.8 9.9±0.3 9.8±0.8 -13.1±0.7 -1.4±1.1 1.9±0.8 -4.8±1.1 0.4±0.9 2.8±0.9 

2 13.6±0.6 16.3±0.6 9.9±0.2 10.6±0.6 -13.0±0.5 -1.8±0.8 1.8±0.6 -4.4±0.8 -0.8±0.7 1.8±0.7 
2a 13.9±0.5 17.7±0.5 10.2±0.1 11.2±0.5 -13.1±0.4 -2.1±0.7 1.5±0.5 -4.0±0.7 -0.5±0.6 1.3±0.6 

3 14.2±0.5 19.3±0.5 10.5±0.1 12.0±0.5 -13.3±0.4 -2.4±0.6 1.4±0.5 -3.5±0.6 -0.5±0.6 0.9±0.6 
3a 15.0±0.5 20.1±0.4 11.2±0.1 11.9±0.5 -13.9±0.4 -3.6±0.6 0.9±0.5 -3.3±0.6 -0.0±0.6 0.7±0.6 

4 17.5±0.5 21.4±0.5 12.3±0.1 11.0±0.6 -13.6±0.4 -4.0±0.6 1.1±0.6 -2.7±0.6 0.3±0.6 1.5±0.6 
4a 20.0±0.5 23.8±0.5 13.3±0.1 10.9±0.6 -13.7±0.4 -3.9±0.6 2.4±0.6 -2.8±0.6 0.6±0.6 1.9±0.6 

5 22.5±0.6 26.3±0.6 14.8±0.2 11.5±0.7 -14.3±0.4 -5.3±0.6 2.2±0.7 -2.6±0.7 0.0±0.7 1.2±0.7 
5a 26.2±0.7 32.9±0.7 17.1±0.2 12.0±0.8 -14.8±0.5 -6.2±0.8 1.5±0.8 -1.4±0.8 -0.6±0.8 1.3±0.8 

6 28.8±0.9 39.4±0.9 19.2±0.2 12.5±1.1 -14.8±0.7 -5.6±1.0 2.2±1.1 -0.7±1.0 -0.8±1.0 2.1±1.0 
6a 33.5±1.2 53.1±1.1 23.8±0.3 13.0±1.5 -14.0±1.0 -8.6±1.3 1.2±1.4 0.0±1.4 1.9±1.4 2.5±1.4 

7 37.2±1.4 61.8±1.3 25.7±0.4 10.9±1.7 -15.5±1.1 -8.3±1.5 4.1±1.6 2.0±1.6 0.7±1.6 1.4±1.7 
7a 34.9±2.1 61.5±1.8 23.2±0.6 10.3±2.6 -16.2±1.7 -4.4±2.2 6.4±2.4 3.9±2.4 -0.7±2.4 -1.6±2.5 

8 24.2±1.6 39.9±1.4 15.1±0.4 16.7±1.8 -13.8±1.2 -5.5±1.8 4.4±1.7 3.6±2.1 0.7±1.8 -2.1±1.9 
8a 19.6±1.2 32.5±1.0 12.5±0.3 14.2±1.3 -15.4±0.9 -7.0±1.3 4.9±1.3 1.5±1.5 0.6±1.3 1.5±1.4 

9 17.5±1.1 31.7±0.9 11.0±0.2 12.8±1.2 -17.3±0.8 -9.6±1.2 4.3±1.2 1.4±1.3 -0.3±1.2 1.7±1.2 
9a 16.1±1.1 29.7±0.9 9.6±0.2 13.8±1.3 -17.1±0.9 -9.8±1.2 4.6±1.2 1.6±1.3 -0.2±1.3 1.3±1.3 

10 12.6±0.9 25.2±0.8 8.4±0.2 12.9±1.1 -16.6±0.8 -9.0±1.1 3.1±1.1 0.5±1.2 0.6±1.1 3.8±1.1 
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Here (Uo, Vo, Wo) is the solar motion with respect to the 
local standard of rest (LSR), which we take from DB98 
to be (10.00 ± 0.36, 5.25 ± 0.62, 7.17 ± 0.38) kms" 1 . In- 
serting (25) into equation (16), we can solve for Wq and 



■ J aip 



in two different ways yielding 



- — — Mot 



J acp 



TTw 



Uo— 



V 



Wo 



J a<p 



= w { 



o: 



(26a) 
Vo. (26b) 



4-2.3. Intrinsic Colors? 

We have mentioned before that extinction is partly to 
blame for the non-uniformity in the distribution of stars. 
In addition, interstellar dust reddens the intrinsic colors 
of stars. To estimate these intrinsic colors, wc assume 
that all stars of a given color in our sample have the 
distance equal to the mean parallax (26b). With the 
standard extinction law (Ry = 3.1) the intrinsic color 
can be approximated as: 



(B-V)o = (B-V) ohs -E B _ v 

1 x d 

« (B-V) ohs -— — , 
riv 



(27) 



where we use the subscript for colors corrected assum- 
ing an average extinction of 1 mag per kpc (e.g., Chen 
et al. 1998). In Table 1, we include these estimates of 
the intrinsic colors, as well as an estimate for the average 
distance (1/Tw). 

Red clump stars with intrinsic B — V G [0.9, 1.1] have 
absolute magnitudes similar to A-type main-sequence 
stars with intrinsic B — V G [-0.1,0.1] (e.g., DB98), so 
that the average distance to the A-type and red-clump 
stars should be similar. Inspecting Table 1, it is obvious 
that the no-extinction case does not conform to this this 
expectation at all. On the other hand, when using the 
extinction-corrected colors, (B — V)o, to determine the 
observed color ranges for the two sub populations, we 
find almost identical distances for the red clump and A 
type stars. We thus conclude that, in the average, the 
extinction corrected colors are close to the intrinsic col- 
ors of our target stars. For the remainder of this paper 
we assume that (B — V)o is a good approximation for 
the intrinsic color. But note that our conclusions do not 
critically depend on this assumption. 

5. Discussion of the Results 

In this section, we discuss only the results obtained 
with K-cr-clipping. These have smaller errors than the 
results obtained using all stars and there are only minor 
systematic deviations between the two sets. 

5.1. Initial Considerations 

Before trying to understand and interpret the results 
derived from the proper motions, we must be aware of 



the kind of stars we are dealing with in the various color 
bins. As discussed in §§2.1 and 2.5, it is well known that 
the kinematics of stars changes systematically with age: 
the age- velocity relation (AVR) . A critical age for a stel- 
lar population is 1.5-2 Gyr. The kinematics of youngest 
stars still carry a significant imprint of the initial condi- 
tions, while older stars have had time to reach an equi- 
librium with the large-scale potential of the Milky Way 
(e.g., Mayor 1974; Gomez & Mennessier 1977). Since 
we are interested primarily in the large-scale proper- 
ties of the Galactic potential, it is important to be able 
to estimate the ages of the stars in our samples. The 
ACT data-base does not allow for sophisticated age esti- 
mates of individual stars, but we can estimate the frac- 
tion of "young" stars in each of the color bins we use. 
Stars bluer than (B — V)o = 0.25 (bins l-5a) have main- 
sequence lifetimes smaller than 1.5 Gyr, and we expect 
these stars to exhibit kinematics appropriate for young 
stellar populations. 

5.1.1. Sample Properties 

One can estimate the fraction of stars younger than 
1.5 Gyr if we assume a rapid post-main-sequence evolu- 
tion and a constant star-formation rate (SFR) over the 
last several Gyr. For stars with a main-sequence lifetime 
of 8.4 Gyr [bin 7, (B-V) ~ 0.6], -18% of the stars will 
be younger than 1.5 Gyr. 

Keeping in mind that the fraction of young stars de- 
creases along the main sequence, we might interpret the 
gradual decrease in non-uniformity of the stellar num- 
ber density (Figure 3) as a decreasing fraction of young 
stars. In Figure 4 we present a measure, 5n(£), of the 
non-uniformity of the number-density distribution. In 
fact, we can use Sn(£) as a proxy for the fraction of 
young stars. 

If we approximate that Sn(£) has just two distinct 
values 5n{t)Y and Sn(£)o for 'young' (t < 1.5 Gyr) and 
for 'old' stars, we can estimate the fraction of young stars 
as 

5n(£) - Sn(£) 



h 



5n(£) Y - 5n{£) ' 



(28) 



For the young stars, we can take the weighted average of 
the bins with (B-V)o < 0.25, yielding Sn(£) Y = 0.126. 
For the old stars, we can use fy = 0.18 for bin #7 and 
invert equation (28) to obtain Sn(£)o — 0.024. We ap- 
ply equation (28) to estimate the fraction of young stars 
among the giants. The results are tabulated in Table 1 
for the case of 1 magnitude extinction per kpc. (If no 
extinction correction is used, the fy values decrease by 
~25% over the values listed.) 

Thus, at the blue end, the stars are both young and 
bright, with ages up to to a few rotation periods of the 
Galaxy, and distances out to 2 kpc. For ever redder 
stars up to B — V ~ 0.7, there is a gradual change to 
fainter, hence nearer, and on average, older stars (cf. 
Table 1). As a consequence, the internal kinematics and 
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the averaging volume of the stars change; the first due 
to the kinematics-age dependence, the latter simply be- 
cause of the distance-color relation on the main-sequence 
(Table 1). 

Beyond B~ V ~ 0.7, there is a more abrupt change 
of stellar properties of the sample, because non-main- 
sequence stars take over to dominate the sample. These 
stars are both brighter and, on average, younger than 
the dwarfs, changing again both the kinematics and the 
averaging volume. If we take the non-uniformity as an 
age indicator, Figure 4 clearly shows that the stars with 
B — V w 1 are younger than both the bluer dwarfs and 
the redder giants. In fact, the 5n{tj value at B — V = 1 
equals that of A- type stars, suggesting an average age of 
only a few hundred million years for stars in this color 
range. 

Using the (B — V)q colors listed in Table 1 and a color 
magnitude diagram representative of the Solar neighbor- 
hood (e.g., DB98, their Fig. 1), we identify the stars in 
color bins 7a to be sub-giants. Bin 8 is a mixture of 
sub-giant, giant and red-clump stars. Red-clump stars 
dominate bins 8a, 9, and 9a. Only the last color bin, 
# 10, predominantly comprise red-giant stars. 

Our Sn(£) analysis above indicates a significant frac- 
tion of young stars in all red color bins (see Table 1). 
This can be understood in the context of ongoing star- 
formation activity in the Solar neighborhood and the 
theory of stellar post-main-sequence evolution (Seidel, 
Demarque & Weinberg 1987; Cole 1998; Girardi et al. 
1998). 

It is worth mentioning that at (B — V)o > 1, the differ- 
ence in luminosity between dwarfs and (sub) giant stars 
is larger than about 3 mag such that the number of red 
dwarfs beyond that color is negligible, in particular after 
K-cr-clipping has been applied. 

5.1.2. Expectations for the Kinematics 

From our considerations in §§2.1 to 2.5, we expect 
the changes in kinematics and averaging volume to be 
reflected in the the Oort constants measurable for the 
stars. Thus, already without the unpleasant effect of 
mode mixing, we expect the Oort constants to change 
gradually blue- ward of B — V ~ 0.7 and red- ward of 
B — V ~ 1, while there might be a more abrupt change 
between these two color ranges. 

An important question will then be: which stars give 
us the "true Oort constants"? The young blue stars 
(t~ms ^ 1-5 Gyr) are not likely to be in equilibrium as 
they still exhibit kinematics associated with their birth 
places. Thus, their streaming velocity field likely de- 
viates systematically from that created by closed or- 
bits, and their distances and velocities are correlated. 
Both effects cause systematic errors when interpreting 
the m = and 2 coefficients as the Oort constants. The 
dwarfs at intermediate colors (0.25 < (B — V)o < 0.7) 
probe a very local volume, the secular parallax is esti- 



mated to be less than 500 pc, and it is likely that their 
measurable streaming field is affected by local anomalies 
(§§2.4 to 2.5). 

The red-clump comprises a mixture of stars of various 
ages, which judged from their longitudinal distribution 
is quite affected by spiral arms and thus subject to the 
same objections as the blue stars. Only samples of giants 
redder than (B — V)o w 1.2 (our color bin #10) might 
be both old enough and distant enough to be unaffected 
by non-equilibrium effects or local anomalies. 

5.2. The Fourier Coefficients Measured 

We plot in Figure 5 as function of intrinsic color 
(B — V)q the results for the Oort constants, the inverse 
mean parallax, the asymmetric drift, A — B, and the 
m = 3,4 terms of Jl t * before (open circles) and after 
(solid squares) mode-mixing correction. 

Obviously, there is a significant difference between the 
raw and mode-mixing corrected values. From our dis- 
cussion in §5.3.2, we expect the "true" values will be 
close to the mode-mixing corrected values, but we can- 
not entirely exclude the raw values as a possibility. We 
will consider the difference between the raw and mode- 
mixing corrected values as an upper limit to the system- 
atic error involved. Regardless of this difficulty, several 
inferences can be made. 

First, there is an obvious discontinuity in kinematic 
properties at (B — V)o ~ 0.75. This discontinuity is re- 
flected in all Fourier coefficients of Jl^ raw or mode- 
mixing corrected. This is unlikely to be caused predom- 
inantly by a difference in distance, since there are no 
significant correlations between the Fourier coefficients 
and 7f _1 (not shown). Presumably more important is 
that these stars span the range in ages where the gra- 
dient in the age-velocity relation is large. The unstable 
nature of the Fourier coefficients in the red-clump region 
illustrates that it is very difficult to determine the true 
value of the Oort constants if the ages of the tracer stars 
are ill-determined. Second, the asymmetric drift and 
A — B are only weakly affected by the systematic errors. 
A — B w 26 for blue stars (uncorrected for mode-mixing), 
which is in nice agreement with A — B = 27 measured 
from Cepheids by Feast & Whitelock (1997). There is a 
trend towards larger values for redder stars. At the red 
end, we again notice the strong discontinuity in the red- 
clump region (B — V)o = 0.75 ± 0.15. The red giant bin 
(#10) has A - B = 26.7 ± 0.7 and 29.5 ± 1.4 for the raw 
and mode-mixing corrected values, respectively. These 
values are consistent with a Galactic circular frequency 
of 28 ± 2 as derived from the proper motion for Sgr A* 
(Reid et al. 1999; Backer & Sramek 1999). 

Third, the values for A before and after mode-mixing 
correction differ by just la for the red giant stars: A is 
between 14.2 ± 0.6 and 12.9 ± 1.1. 

Fourth, except for the very blue stars, there is a large 
systematic uncertainty for B. The red giant bin has B 
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Fig. 5. — Inverse mean parallax, the asymmetric drift from equation (26b), Oort's constants (left) and A — B, as well as the higher- 
order terms in the Fourier expansion of ~p^* (right) obtained for the 19 color bins defined in Table 1, after K-cr-clipping, from top to 
bottom. Both the raw results (circles) and the mode-mixing corrected values (squares) are presented. The vertical dotted lines indicates 
(B-V) = 0.75. 



values between -16.6 and -12.5. 

Fifth, there are several clear evidences for deviations 
from axisymmetric equilibrium: non-zero C and sin 31 
terms (the latter only for blue stars). Note, in particular, 
that after mode-mixing correction C has the same sign 
for all stars. 

Finally, while the asymmetric drift is increasing with 
color independent of whether it has been derived from 
the raw or mode-mixing corrected coefficients, the curves 
derived from the U and the W motion differ clearly (not 
shown). In the latter case, the increase is more gentle 
and reaches only v aip m 16 krus , whereas the [/-derived 
coefficients yield v aLp 2-4 km s _1 larger. Note that we do 
not expect an exact correspondence with the asymmet- 
ric drift derived by DB98 because our samples comprise 
mixes of sub-populations that differ from the local Hip- 



parcos sample, in particular at the red end. 

5.3. Mode Mixing 

5.3.1. Evidence for Mode Mixing 

In Figure 6, we plot the Fourier coefficients (cv, s 7r ) TO , 
derived from the vertical proper motions utilizing equa- 
tion (18), versus the inverse mean parallax estimated 
from the average vertical proper motion via equa- 
tion (26b). Clearly, these coefficients deviate signifi- 
cantly from zero proving that mode mixing is present 
and non-negligible, i.e. that the Oort constants ob- 
tained without correction are systematically in error. 
There is a clear dichotomy between main-sequence stars 
(B — V < 0.8; open symbols) and non-main-sequence 
stars (B — V > 0.8; closed symbols): most mode-mixing 
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Fig. 6. — Mode-mixing coefficients (defined in equation (14)) as 
estimated from equation (18) for the bins 1 to 7 (open) and 7a to 
10 (solid) plotted against the inverse mean parallax estimated via 
equation (26b). 



coefficients are similar within each color group but differ 
between them. This indicates different spatial distribu- 
tions, necessitating separate analyses. 

5.3.2. Is our Mode- Mixing Correction Correct? 

We have tried to correct for the mode mixing by em- 
ploying the vertical proper motion and the technique de- 
scribed in §3.2.1 and §3.3.2. However, as discussed in 
§2.7.3, there are possible caveats in this method, in par- 
ticular neglecting contributions from higher-order terms 
(cf. §§2.4 to 2.5). One possible check on the consistency 
of the results obtained is a comparison of the secular 
parallaxes estimated from the solar radial and vertical 
motions via equations (26). Figure 7 compares the esti- 
mate Ify 1 , which is affected by mode mixing, with 7r^, 
which is not, before (circles) and after (squares) mode- 
mixing correction has been applied. Obviously, in both 
cases a systematic deviation of 5-10% exists between the 
two estimates, but with opposite signs. An exception are 



"cs 




Fig. 7. — The mean distance estimated from the vertical proper 
motions (via equation (26b)) versus the distance from the radial 
proper motions (via equation (26a)). The raw results are repre- 
sented as circles, the mode-mixing corrected results as squares. 
The size of the symbols increases with (B — V)q. 



the very early-type stars (B — V < 0.2), which are known 
to deviate from equilibrium (DB98). 

Note that while the deviations from the dotted line 
increase with distance, they do so less for the mode- 
mixing corrected results. This indicates that the neglect 
of higher-order terms, which should introduce system- 
atic errors at large distances, cannot have introduced 
significant errors. 

There are two possible explanations for those differ- 
ences. First, the ratio Uq/Wq may be different for this 
sample than for the more local sample of Hipparcos stars 
utilized by DB98. In this case, we cannot make any 
statement, whether our mode-mixing correction works 
or not. 

Second, if Uq/Wq is equal the DB98 value, we can 
determine an alternative mode-mixing solution in which 
we force equality between nu and ttw ■ In this case, the 
radial and tangential proper motions are given by 

TT- U ° ~ 



s 7r2 x HuO.W ~ C H 
1 + C n2 



from equations (26) and (15a), respectively (we have 
dropped the latitude dependence). When using these re- 
lation for ~p uo and ~Pvo an d solve the mode-mixing equa- 
tions (15a) directly, we find that the so-determined Oort 
constants do not differ substantially from our previous 
mode-mixing corrected values. This indicates that the 
slight difference between du and dw does not signify a 
substantial problem for the mode-mixing scenario. 
Since there is overwhelming evidence for the reality 
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Table 5 

Properties of some critical sub-populations 



(B-V)o d f Y T B r B m 



< o 


0.8 


15 


6 


100 0.1 


-12 


-13 


0.6 


0.3 


36 


15 


30 8 


-10 


-16 


1.2 


0.8 


38 


17 


30 8 


-12 


-16 



Note. — We list approximate values for the intrinsic color, dis- 
tance, velocity dispersion, asymmetric drift, percentage of young- 
stars, approximate average age [Gyr], and the raw and mode- 
mixing corrected Oort's B (B r and B m , respectively). 

of the mode-mixing effect (Figure 6) and the details of 
the mode- mixing procedure appear to be irrelevant, we 
conclude that the mode- mixing solutions arc robust. In 
particular, we are confident (and hope to have convinced 
the reader, too) that the remaining systematic errors 
of the corrected results are signigicantly smaller than 
for the raw values. In our subsequent analysis, we will 
concentrate on the mode-mixing corrected values, but 
will also show the raw results for comparison. 

5.4. Can we make sense of these results? 

In §2, we discussed various potential origins for de- 
viations of the measured Oort constants from their 
"true" values. These deviations, which can be several 
kms _1 kpc _1 , originate from departure of the Milky 
Way from a smooth axisymmetric equilibrium, and may 
depend on velocity dispersion and mean depth of the 
stellar population considered. Variations of this order 
are seen in our results in Figure 5. While a detailled in- 
terpretation of the measured variations of the Oort con- 
stants is beyond the scope of this paper, we may non- 
theless examine whether we can single out a dominant 
cause 8 . 

To disentangle the factors that may contribute to the 
variation of B with v aip , we first compare three sub- 
populations in Table 5. The bluest and reddest stars 
have similar distances but very different ages, and hence 
also different velocity dispersions. As the third, we take 
the reddest main-sequence stars at (B — V) ~ 0.6, which 
are much nearer. 

In fact this is the only relevant property in which they 
differ significantly from the red giants, see Table 5 and 
Figure 5, in particular if one considers the mode- mixing 
corrected values. A straightforward interpretation is 
that differences in sample depth are unimportant in caus- 
ing the differences in the observed Oort constants. This 



We like to mention at this point that any systematic error in the 
cataloged proper motions also adds to these variations. However, 
the variations persist when performing the same analysis on the 
Tycho Reference Catalogue (Hog et al. 1998) or the Tycho-2 Cat- 
alog (Hog et al. 2000b). This latter catalog is based on the same 
data as the ACT, but its proper motions are derived in a differ- 
ent manner, so that it seems unlikely that the proper motions are 
systematically in error (see also Urban, Wycoff, & Makarov 2000). 



in turn implies that small-scale wiggles of the Galactic 
velocity field cannot be important either, for otherwise 
we would expect significantly different results for red MS 
stars and giants. 

We can use the asymmetric drift as a proxy for the 
mean age and velocity dispersion of a stellar population 9 . 
Similarly, we may use the overabundance Sn(£) of stars 
in directions of spiral-arm tangents as a proxy for the 
fraction of young stars. In Figure 8, we plot the Oort 
constants measured, including the m = 3 Fourier coef- 
ficients, against our estimate (26b) for the asymmetric 
drift (left) and 5n(£) (right). The raw data (circles) show 
only marginal trends with v av , and perhaps even some 
dichotomy between red main-sequence (small symbols) 
and giant populations (large symbols). On the other 
hand, the mode-mixing solutions do not show such a di- 
chotomy but rather tight linear relations between the 
Oort constants and the asymmetric drift. 

A dichotomy is clearly present in the plots versus 
5n(£), which arises because the red-clump stars at 
8n(£) ~ 0.13 have values similar to those of the red 
giants and unlike those of the blue MS stars at the same 
5n(£). Remarkable, however, is that the red MS stars 
and red giants at 8n(£) k, 0.04 show rather similar val- 
ues, despite that fact that their average distances differ 
by a factor > 2. This is true even for the high-order 
terms, which argues for their non-zero values not being 
created by sample-depth effects. 

In §2.1 we saw that there should be almost no corre- 
lation between asymmetric drift and Oort's B, for the 
case of axisymmetric equilibrium (and constant sample 
depth). Thus, the very presence of such a correlation 
with v a , v as well as non-zero C and s^3, argue strongly 
for non-axisymmetry to be the dominant origin of the 
observed differences of the Oort constants between stel- 
lar subsamples. 

Alternatively, non-equilibrium effects may play a role. 
However, one would then expect a somewhat erratic be- 
haviour, instead of the clean trends seen in the left plot 
of Figure 8 and also in Figure 5 for the main-sequence 
stars (i.e. from (B-V) = -0.3 to 0.6). 

To summarize, the variations of the Oort constants 
between the subsamples originate most likely in de- 
viations from axisymmctry, while non-equilibrium ef- 
fects and small-scale wiggles (deviation from smooth- 
ness) seem less important (as do any effects that rely on 
variations in sample depth). This in turn implies that 
any interpretation of the Oort constants in terms of the 
properties of the underlying Galactic potential, the very 
motivation to undertake studies like this, see §1, cannot 
be as simple as in Oort's days. 



Recall that stars obey an age-velocity dispersion relation (<r% oc t; 
e.g., Jenkins 1992) as well as Strombcrg's asymmetric drift relation 
{al oc v alp , e.g. DB98), i.e. v av oc t oc a 2 t . 



19 




10 15 

v av ,w [ km s_1 ] 



20 

15 

10 

5 
-5 
10 
15 



'cq 30 

I 

^ 25 

20 
1 

5 

*o 

-5 
-10 

5 

"< 

-5 

5 

-5 



l Co 



h-|H 



I 1 1 h 




n — xi — 1 — r 




0.1 



0.2 



0.3 



Fig. 8. — Plotted as a function of the asymmetric drift (left) and the overabundance Sn(£) of stars in direction of spiral-arm tangents 
(right), we present A, B, A — B, C, C13, and 5^3, from top to bottom. All "Oort constants" depend clearly upon v a A, most significantly 
so in the mode-mixing case (squares). The symbol size grows with increasing (B — V)q 



Summary and Conclusion 



the Oort constants are denned as the divergence (K), 
vorticity (B), azimuthal (A) and radial (C) shear of the 
local stellar streaming field of the Milky Way in the limit 
of vanishing random motions, where all stars move on 
closed orbits, i.e. circular orbits in case of axisymmetry. 
The importance of the Oort constants derives from the 
fact that the dynamics of these closed orbits is directly 
related to the Galactic gravitational potential, a relation 
that becomes particularly simple in the axisymmetric 
case. 

The longitudinal proper motion of a star with parallax 
7r and radial and azimuthal velocity with respect to the 



Sun, U and V , may be written as 

Ji r = -Rc{tt{V + iU)e lt }. 



(29) 



Assuming that the stellar kinematics and parallaxes are 
uncorrelated, one finds for the mean longitudinal proper 
motion 

■p tt = -Re {if (V + iU)e l1 }. (30) 

The spatial variations of U and V are given by the Oort 
constants, and lead to a |m| = 1 harmonics % %ml in their 
Fourier expansion. Together with the e tl in (30) this re- 
sults in m — 0, 2 harmonics with amplitudes given by 
the Oort constants A, B, and C . A similar harmonic de- 
pendence is exhibited by the mean radial velocity times 
parallax, -kv t . However, stellar parallaxes and radial ve- 
locities are difficult to measure and thus the Oort con- 
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stants are most commonly determined from their effect 
on the stellar proper motions. 

6.1. Mode Mixing and Other Problems 

The effect of the Oort constants on the stellar proper 
motions is comparably small. For nearby stars (within 
~lkpc), it is much smaller than the contributions from 
random stellar motions (i.e. dispersion of U and V in 
equation (29) and the reflex of the solar motion (i.e. the 
lowest order in equation (30). Therefore, large proper 
motion surveys are necessary to extract the Oort con- 
stants with reasonable accuracy. There are various, 
mostly known but neglected, problems arising in this 
procedure, which are summarized in the footnotes of Ta- 
ble 6. 

A fundamental, hitherto apparently overlooked, prob- 
lem in measuring the Oort constants from proper motion 
data is what we call mode mixing (point vii in Table 6). 
In fact, the cause of the problem is rather similar to 
the very effect one is after. As a spatial variation of 
(U, V), described by the Oort constants, contributes to 
the m = 0, 2 harmonics in the Fourier expansion of /J^* , 
so does a variation of n. This contribution is indistin- 
guishable from that due to the Oort constants itself. 

For a typical situation, already a variation in the mean 
parallax of only 10% results in a contribution of a few 
kms _1 kpc _ , larger than any other source of uncer- 
tainty in the Oort constants. Such a variation of W is not 
anticipated from a smooth exponential disk. However, 
that seems to be a bad description of the actual situation 
for a typical stellar sample. First, the stellar distribu- 
tion is inhomogeneous, in particular for early-type stars, 
which are predominantly situated in spiral arms. Sec- 
ondly, even if the underlying density is rather smooth, 
extinction inevitable leads to significant inhomogeneities 
in any actual stellar sample. 

In order to correct for the mode mixing, one needs 
an unbiased estimate of the variations of W with I. For 
low-latitude samples, such an estimate may be provided 
from the latitudinal proper motion m, the mean of which 
at b = should be given by JZ b — n W. Here, W is the 
vertical stellar motion with respect to the Sun. If (1) ir 
and W are uncorrelated, and (2) W is constant with £, 
then we can directly measure the relative variations of n. 
A local stellar warp would invalidate both these assump- 
tions, and, hence, one would best not rely on them, but 
use an independent and unbiased, though not necessarily 
very accurate, distance estimate for the individual stars. 
Unfortunately, such an estimate is not currently avail- 
able - note that the photometric parallax is ill-suited for 
this purpose. 

There is also a more technical issue that we want to 
emphasize. The standard way to determine the Oort 
constant from proper motion surveys used to be a least- 
square fit of the data to a functional form that just ac- 
counts for the reflex of the solar motion and the Oort 



constants, i.e. does not allow for higher-order Fourier 
modes. This essentially is a parametric fit, and as such 
will give biased estimates for the parameters if the data 
are not well described by the functional form fitted. In 
our analysis, we allow for high-order terms, and found 
indeed that these are non-negligible. We also restricted 
our analysis to low latitude stars for three reasons: (i) 
restricting the latitude to a narrow range, we do not 
need to account for possible latitudinal variations of the 
mean proper motions; (ii) only at low b can we use the 
latitudinal proper motion to correct for mode mixing; 
and (iii) the Oort constants are strictly defined only for 
closed, i.e. planar, orbits. Moreover, most disk stars in 
sufficiently deep samples are predominantly at low lat- 
itudes anyway. A disadvantage of low latitude stars is 
that extinction is largest and most patchy in this region, 
so that the variation of W with longitude is likely to be 
maximal here. 

6.2. Longitudinal Number-Density Variations 

Before analyzing the proper motions of the low- 
latitude (|6| < 5.73°) stars from the ACT/Tycho-2 cat- 
alogs, we first considered their longitudinal distribution 
as a function of color. For the bluest stars in the sample 
(B — V w —0.1), this distribution shows three narrow 
peaks in the direction of the nearby spiral-arm tangents. 

Next, we quantified the over-abundance of stars at 
the longitudes associated with these three peaks. This 
over-abundance may be interpreted as a measure for the 
importance of stars in spiral arms. In fact, we use the 
relative over-abundance Sn(£) to estimate the fraction of 
young stars in our color bins. These estimates compare 
well with estimates based on the star-formation history. 
With increasing color along the main sequence, 5n(£) 
decreases, but never diminishes completely. 

At colors redder than (B — V)o ~ 0.75, sub-giant, red- 
clump, giant and super-giant stars dominate the sample. 
Interestingly, at about (B — V)o = 0.9, i.e. for red-clump 
stars, the over-abundance of stars at the peak-longitudes 
increases back to the value for early-type stars as blue 
as (B — V)o w 0.25. This can be explained by the about 
equal duration of the horizontal-branch phase for stars of 
all masses <2M Q , making the red clump a strange mix- 
ture of stars of all ages. As a consequence, any external 
observer of the Milky Way using some red pass band 
which allows for red-clump stars to contribute signifi- 
cantly, would see a spiral pattern considerably stronger 
than that of the underlying old-disk population, which 
essentially represents the stellar mass distribution. 

For yet redder colors, i.e. beyond (B — V)o ~ 1.2, 
our measure for the importance of spiral structure 
drops again, but not quite as low as for the dwarfs 
at (B — V)o ~ 0.6. This is because the duration of 
the giant-branch phase is shorter for more massive and 
hence younger stars, such that the giants are rather old 
on average. 
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Table 6 

The Oort Constants: true and measurable 







symbols 


definition and description 


(1) 


A, 


B, C, K 


the Oort constants: 

divergence, vorticity, and shear of the (hypothetical) velocity field v(x) due to the 
closed orbits supported by the Galactic potential 


(2) a 


A 


B, C, K 


the best we can hope to get: 

divergence, vorticity, and shear of the streaming velocity field v(x) of a group of stars 


(3) b 


A 


B, c, k 


what we can actually measure: 

Fourier coefficients of the proper motions measured for a group of stars 



a Possible reasons for differences between (1) and (2): 

(i) For young stars: moving groups and other non-equilibrium effects like spiral arms lead to unpre- 
dictable deviations of v from a closed-orbit streaming field. 

(ii) Local anomalies in the streaming field are reflected in the Oort constants if the sampling volume 
is too small. This is mainly affects sub-populations with low velocity dispersion making them 
susceptible to small-scale variations in the Galactic force field. 

(iii) For old stars: v deviates from closed-orbit streaming by the asymmetric drift. For the axisymmetric 
case, this effect can be estimated using Sternberg's asymmetric-drift relations: it reduces A by up 
to 3kms _1 kpc _1 , but hardly changes B. 

b Possible reasons for differences between (2) and (3): 

(iv) Correlations between the stellar parallaxes and velocities, which may occur for stars associated 
with spiral arms or a local warp, invalidate the basic assumption underlying the Fourier approach. 

(v) Terms of higher order than linear (= the Oort constants) in the Taylor expansion of the streaming 
field become increasingly important with ever deeper samples. 

(vi) Discontinuities of the streaming field at the OLR of the Galactic bar render the Oort constants ill- 
defined. If the sampling volume contains such places of resonance, the volume-averaged streaming 
field deviates systematically from the local field. 

(vii) Mode mixing: variations of W with t in conjunction with the solar reflex motion lead to contri- 
butions to the proper motion that are indistinguishable from the Oort constants and up to a few 
kms -1 kpc -1 in size. 
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6.3. The Oort Constants 

The significant non-uniformity of the longitudinal 
stellar distribution strongly hints towards inhomogene- 
ity in the spatial distribution of the sample stars. Thus, 
longitudinal parallax variations and hence a confusion of 
the Oort constants by mode mixing with the solar reflex 
motion is expected and observed. In the previous sec- 
tions we have presented evidence that the mode-mixing 
corrected values are likely to be close to the "true" Oort 
constants. We thus present the mode-mixing results be- 
low. 

The absolute values of the Oort constants are small 
for young, blue stars. Extrapolating to zero asymmetric 
drift we find (in kms -1 kpc -1 ). 



young MS stars: 



A = 


9.6, 


B = 


-11.6, 


A-B = 


21.1, 


C = 


0.4, 


S£,3 = 


-6.7, 


Q,3 = 


-0.1, 



(31) 



with internal errors of about 0.5 km s -1 kpc - . The old- 
est stars in our samples, giants with (B — V)o > 1-2, have 
significantly larger Oort constants. At v aip — 18 we get: 



old red giants: 



A = 
B = 
A-B = 
C = 

se,3 = 
ct.z = 



15.9 
-16.9 
32.8 
-9.8 
2.6 
4.3 



(32) 



with internal errors of about 1.2kms _1 kpc -1 , and 
where 2 kms -1 kpc -1 has already been added to the A 
value of the giants in order to correct for the reduction 
by the asymmetric drift (cf. §2.1). The external errors 
are harder to estimate 10 . 

Our values for A and B for the old stars are some- 
what unusual, which is largely due to our mode-mixing 
correction. Using the raw Fourier coefficients whereby 
ignoring this effect, one finds for the giants more conven- 
tional numbers, in particular A raw ss 14, £? raw ~ —12.5 
and C raw « 0. 

The values for the Oort constants for the young and 
old stars derived with the mode-mixing corrections are 
not inconsistent with previous determinations. In part 
this arises due to the fact that most previous work has 
applied parametric methods to derive the Oort con- 
stants, in part because published results still have sig- 
nificant uncertainties, and partly because all these val- 
ues are uncertain due to the various systematic effects 
discussed above. 



10 If we assume that the difference between the raw and mode-mixing 
corrected values are indicative of such external errors, we find: 
[8A,Ss e:i ,SB, SA-B,Sc ii3 ,SC] » [0.1,0.5,4.4,4.3,4.3,9.8] fs [0.6, 
19, 26, 13, 100, 100]% for old stars, and almost zero for young 
stars. 



A relatively large value for f2 has been found from the 
longitudinal proper motion of Sgr A* (if Sgr A* is at rest 
in the Galactic center). In that case, £1 = A — _B = 28±2 
(Backer & Sramck 1999; Reid ct al. 1999). Even larger 
values for A — £?(« 31.5) have been reported by Mendcz 
et al. (1999) and Miyamoto & Zhu (1998), based on ab- 
solute proper motions derived from the Southern Proper 
Motion Program and the Hipparcos Catalogue, respec- 
tively (note, however, hat these two works have employed 
the classical parametric /i£ = B + A cos 2£ description of 
Galactic rotation). All of these results are consistent 
with our determinations for old stars (32). 

6.4. What can be learned from the Oort con- 
stants? 

To our knowledge, the linear dependence of the Oort 
constants on asymmetric drift is a new discovery. It is 
worth noting here, that while such a dependence was in 
principal expected already by Oort, our measurements 
are inconsistent with expectations from axisymmetric 
models for the Milky Way. 

The non-zero §1,3 value for young stars clearly shows 
that non-axisymmetry is present amongst these stars, 
which given their affiliation with spiral arms is not sur- 
prising at all. While for older stars % 3 diminishes, an- 
other indicator of non-axisymmetry, Oort's C, becomes 
very significantly non-zero. For other sitations than ax- 
isymmetry, predictions for the values expected for the 
Oort constants are considerably involved and require de- 
tained modelling. This is true cspcccially, if effects of the 
sample depth and velocity distribution are to be taken 
into account. 

As indicated by the imprints of the Galactic bar on 
the local velocity distribution measured from Hipparcos 
data (Dehnen 1998), the Sun lies just outside the outer 
Lindblad resonance (OLR) of the bar (Dehnen 2000). 
This is exactly the place where the bar is expected to 
most strongly distort the local kinematics. In particular, 
stars from inside the OLR should produce C < 0, and in 
a deep sample the number of such stars should be larger 
than in a very local volume. Indeed, when analyzing 
only the stars brighter than Vt — 10.5, C rises to —7.5, 
i.e. the effect is weaker. Thus, the value of C might 
arise from contributions of stars from inside the OLR. 
Miihlbauer & Dehnen (2003) have analysed the effect of 
the Galactic bar on the stellar kinematics in the outer 
disk and find that larger a velocity dispersion, as for 
the red giants, tends to shift the "effective radius" of 
resonance outwards. They also find weaker distortions 
to the local C than our result. 

Another intriguing possibility is that asymmetries in 
the stellar distribution function arising from a (or more) 
satellite merger, like that of the Sgr dwarf, contributes 
to the apparent anomalies. Note that there appear to be 
asymmetries in the local number of thick disk/inner halo 
stars (Larsen & Humphreys 1996; Parker, Humphreys, 
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& Larsen 2001). 

Further work is needed to investigate whether such 
values can be explained when (1) considering sample- 
depth effects, i.e. simulating the actual data-taking and 
measurement process, (2) residual non-equilibrium ef- 
fects from minor mergers, (3) and accounting for the 
additional non-axisymmctric forcing of triaxial halo, el- 
liptical disk (e.g., Kuijken & Tremaine 1994), or spiral 
arms. This latter option is interesting as the Sun resides 
also close to the co-rotation radius of the spiral pattern 
(e.g., Lin, Yuan & Shu 1969; Mishurov Pavlovskaia & 
Suchkov 1979), where one would expect its forcing to be 
most efficient. 

6.5. Resume 

The very concept of the Oort constants is originally 
based on the idealization of the Milky Way being ax- 
isymmetric with almost circular stellar orbits. Similarly, 
in determining the Oort constants from stellar proper 
motions one has usually relied, at least implicitly, on 
the smoothness of the spatial distribution of the sampled 
stars. At Oort's time, the deviations of reality from these 
ideals have been of second order compared to the internal 
uncertainties due to the low accuracy and sparse sam- 
ples. However, with ever increasing accuracy and larger 
samples, these deviations inevitably become significant. 
Nonetheless, they have been largely ignored so far, re- 
sulting in significant systematic errors. The presence of 
such errors is reflected in the previous determinations of 
the Oort constants listed by Kerr & Lynden-Bell (1986), 
whose scatter exceeds the typical internal errors. 

In this paper, we (1) tried to explain the possible de- 
viations from the above ideal and the resulting system- 
atics, in particular when using large proper motion sur- 
veys, and (2) tried as far as possible both to avoid and 
correct for them in our application to the ACT/Tycho-2 
catalogs. Avoiding uncontrollable systematics lead us to 
exclude the results from either young, high-latitude, or 
nearby stars. Thus, the most reliable results reported 
in the present study are those for the (mainly old) gi- 
ants in the Galactic plane, for which there are two main 
sources of systematics: the asymmetric drift and mode 
mixing. While the first problem was already known to 
Oort (1928), the second is not mentioned in previous 
determinations of the Oort constants from proper mo- 
tions. Our correction for this latter effect using the stars' 
vertical proper motions changes the Oort A, B, and C 
from more conventional numbers to the values reported 
in equation (32) for old stars. 

We hope that possible future space missions such as 
FAME, DIVA, SIM and GAIA 11 will dramatically im- 
prove both the quality and amount of the data, yield- 



FAME:http://www. usno.navy.mil/FAME/ 
DIVA: http://www.aip.de/groups/DIVA/ 
SIM: http://sim.jpl.nasa.gov/ 
GAIA: http://astro.estec.esa.nl/GAIA/gaia.html 



ing highly accurate five- or even six-dimensional phase- 
space coordinates for millions of stars throughout the 
Milky Way. With such data, many of the hitherto un- 
pleasant deviations from Oort's idealized Galaxy will be- 
come clearly apparent and must be dealt with properly. 
Very likely, the classical deductive way of 'measuring' the 
Oort constants directly from the data, becomes hardly 
viable then. Instead, a detailed analysis is needed which 
compares the data with predictions from sophisticated 
dynamical models incorporating the Galactic bar, spiral 
structure, non-equilibrium effects such as moving groups, 
and other deviations from Oort's ideal Milky Way. In 
fact, these deviations, which plague our current determi- 
nation of the Galactic potential, are likely to constrain 
the models significantly leading eventually to a much en- 
hanced understanding of the structure and formation of 
our host galaxy. 

One datum that is very important to the study of 
the dynamics of the Milky Way will not be collected by 
the future space-based astrometric missions. The lack 
of approximate ages of the individual 12 tracer stars will 
seriously compromise the study of the dynamical evo- 
lution of the Galaxy. We have already noted that the 
interpretation of the proper motions of certain kinds of 
stars (the red-clump region) is hampered due our in- 
ability to associate these stars with either young or old 
stellar populations. 

We anticipate that age-related effects will become an 
important limitation to the interpretation of the space- 
based astrometric data. Age information would allow us 
to study many topics related to the formation and evolu- 
tion of the Milky Way, such as (1) the (relative) ages of 
various stellar components, (2) the temporal variation of 
the star formation activity, (3) the minor-merger history 
(through distinct features in the age- velocity relation). 
Such studies would turn the Milky Way into an impor- 
tant benchmark galaxy, not only at the present, but also 
at cosmologically interesting epochs. 

Substantial parts of this work were performed at Ox- 
ford, Heildeberg and Potsdam (WD), and Southampton, 
Rutgers and USNO (RPO). 
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